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ABSTRACT 

Although parts of this paper may be applicable to various 
other location finding systems/ radio direction finding was 
selected for development because of the extensive literature 
available on the subject. Using a modification of the navi- 
gation coordinate system, equations in spherical trigonometry 
are developed by analogy with graphical solution procedures. 
Ellipticity of the earth is corrected for by use of geocen- 
tric latitude in all calculations. The least squares solu- 
tion is found by a combination of second-order gradient 
searches and Golden Section Searches on the corrected sphere. 
A Chi-Square contour confidence region is presented with 
defense against the empty region arguments and with verifi- 
cation of its reliability of target containment. Sensitivity 
of location error to network layout, number of measurements 
available and intercept geometry factors is examined by 
computer simulation, demonstrating the predictive value of 
these statistics. It is shown that the correlation between 
location error and the size of a reliable confidence region 
is too minimal to have predictive value. 
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I. 



INTRODUCTION 



Various procedures have been described for estimating an 
emitter location from a set of measured signal azimuths 
reported to a control station by dispersed direction finding 
(DF) sites. Such procedures range upward in complexity from 
visual estimates off a plot of great circle arcs on a 
gnomonic projection of the operating area of potential 
targets. Other procedures exist for describing confidence 
criteria about the location estimate. The present paper was 
motivated by an interest in confidence regions and in loca- 
tion error analysis. Such studies require either that the 
procedure used be free of approximating assumptions, or that 
the variability of the approximations be included in the 
confidence criterion. The former course was preferred for 
simplicity and consistency of estimates. This paper presents 
an emitter location estimate procedure based only upon known 
DF site coordinates and signal azimuths. The procedure 
described in sections II and III converges quickly to a 
least squares solution and does not require an initialization 
point estimate of the location from external sources (Pope, 
1971) . The procedure does not require a rejection of out- 
liers (McCalla, 1970) in order to converge to the least 
squares solution. The approach is free from assumptions in 
the categories that the earth is flat in the solution region, 
that the earth is a perfect sphere, that great circle arcs 
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are well enough approximated by straight lines of position, 
and that azimuth errors are well enough approximated by 
parallel displacements of such lines of position (Daniels, 
1951 ) . 

Section IV describes a Chi-Square contour confidence 
region procedure. Factors from intercept geometry are 
discussed in section V. Section VI is a description of 
procedures used to simulate the problem for the purpose of 
the location error analyses described in sections VII and 
VIII. An alternative to confidence regions for predicting 
location accuracy is discussed in the concluding section. 
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II . 



MATHEMATICAL DESCRIPTION OF SPATIAL RELATIONSHIPS 



Locations and spatial relationships of interest in direc- 
tion finding and target fixing can be expressed completely 
in spherical trigonometry. The advantage of such expression 
is that each mathematical step has its simple graphical 
counterpart. This paper contains notational and descriptive 
departures from other direction finding models based on spher- 
ical trigonometry. Such notational changes follow those used 
in geodesy for describing locations and spatial relationships 
and those used in textbooks on spherical trigonometry for 
developing equations to express such relationships. Equations 
are restricted to sine and cosine functions only in order to 
simplify calculating first and second order derivatives of 
somewhat lengthy equations. 

A. GLOBAL COORDINATE SYSTEM 

The coordinate system is a slight modification of the 
system used in mapping. Longitude, denoted X, is measured 
from zero at the prime meridian eastward through 360 degrees. 
Latitude, denoted <j> , is measured from zero at the equator 
northward through ninety degrees and southward through minus 
ninety degrees. Azimuth, denoted 0, is measured from zero 
at true north clockwise through 360 degrees. Such a coor- 
dinate system is sufficient to represent any point on the 
globe as well as the great circles which represent ideal wave 
propagation paths. Mathematical descriptions of spatial 
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relationships must account for the longitude discontinuity 
at the prime meridian and for the azimuth discontinuity at 
true north. A great circle arc with azimuth angle indicated 
is sketched on the coordinate system in Figure 1. 




Global Coordinate System 
FIGURE 1 



B. GEOCENTRIC LATITUDE 

Geographical latitude, used in navigation and mapping, 
is based upon astronomical observations and mathematical 
deductions that do not correct for the ellipticity of the 
earth. Consequently, geographical latitude is slightly 
greater in magnitude than the arc on a meridian from the 
equator to a location with vertex angle at the center of the 
earth. Rules of spherical trigonometry can be more accurately 
applied to describe spatial relationships following the con- 
version from geographical latitude to geocentric, or earth- 
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centered, latitude. Assuming that waves propagate along 
great circle arcs, then no information is lost since such 
paths are geometrically independent of radii to the earth 
center from both the signal emitter and from the DF site. 

The maximum discrepancy between the two types of latitude 
occurs at ±45 degrees latitude and is about equal to 11.5 
minutes of arc. Although the conversion equation appears 
in books on geodesy, a derivation could not be found, so the 
equation was derived for reference purposes. 

In Figure 2, a and b represent the semi-major and semi- 
minor axes of the ellipsoid. The point P(x^,y^) represents 
any point on the ellipsoid. The line OP connects P to the 
origin at 0, the earth center. SP is perpendicular to the 



y 




Relationship Between Geographical and Geocentric Latitude 

FIGURE 2 
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equatorial plane, the x-axis. NP is the normal line to the 
ellipsoid at P, and the angle PNS is equivalent to geographi- 
cal latitude. Angle PON is equivalent to geocentric latitude, 
denoted ip. Ellipticity of the quadrant is exaggerated in 
order to distinguish the two angles. 

From Figure 2, it can be seen that 



SP = OS tan ip, and that SP = NS tancf> , 



so that 



. . NS , , 

tanif> = Qg- tan(J) 



( 2 - 1 ) 



From the equation of the ellipse. 



b . 2 2 . 2 

*1 = a <a ' x l> 



(2-2) 



The equation of the normal line NP is given by 



y - y i = 



f T Tx7T (x " x i ) 



(2-3) 



Differentiating equation (2-2) to solve f'(x^), 

- 1 - 

, . b ,2 2, 2 

f (x^ = _ ^ X 1 ( a _ x ;l) 

Substituting in equation (2-3) and rearranging terms. 



y = bkr (a2 - x i> 5 (x - x i> 



2 

a Y 1 

= — 2 (x - x^) , from equation (2-2) 

b x^ 



Then at y = 0, y^ can be cancelled from both sides, and 
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X 



- X 



1 



X 



1 




2 

a 



From Figure 2, x^ = OS, and at y = O, x = ON, so that 
x^ - x = NS. 

NS 

OS 2 
a 



Substituting this last result in equation (2-1) , 



tantjj 



^2 tan<f> 
a z 



(2-4) 



The conversion is usually stated as a function of earth 
eccentricity. 



e = 



1 

a 





Then 



2 

e = 




b 



2 



Therefore , 





tantjj = 



2 

- e , and it follows that 
(1 - e 2 ) tan4> 



(2-5) 



Using e = 0.081813 from the International Ellipsoid of 
Reference (reference 3) , 



tamf> = 0.99327733 tan* 



( 2 - 6 ) 



C. THE POLAR TRIANGLE 



The mathematical description of spatial relationships is 
based upon a polar triangle constructed from the great circle 
arcs (measured by the angle subtended by the arc at the 
earth's center) connecting a DF site, any other point and the 
north pole. Such a triangle is sketched in Figure 3. The DF 
site is located at vertex B. Vertex C is at any other point. 
Vertex A is at the north pole. The arc c is the complement 
of the geocentric latitude of the DF site. Similarly, the 
arc b is the complement of the geocentric latitude of -the 
other point. Angle A is equal to the difference in longitude 
between the DF site and the point. The dotted arc included 
on Figure 3 corresponds to a signal azimuth with angle indi- 
cated by 6 . 



A 




Polar Triangle with Signal Azimuth 
FIGURE 3 
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Angle A corresponds to the smaller of two possible polar 
triangles connecting any two points and the north pole. 
Visually, there is no problem in deciding which polar tri- 
angle to choose. However, when this triangle contains the 
prime meridian, the discontinuity at longitude 360 = longi- 
tude 0 results in assigning the value of the vertex of the 
larger polar triangle when the longitude difference is taken. 
The mathematical model can be forced to correspond to the 
spatial relationships by adding the proviso that if the 
difference in degrees of longitude is greater than 180 
degrees, then let A equal 360 degrees minus the difference. 
This restricts the size of angle A to one hemisphere. The 
absolute value of the difference can be taken in case the DF 
site is east of (greater longitude than) the point. These 
relationships are summarized below with subscript s desig- 
nating the DF site and subscript t designating any test 
point. 



c = 90 - ip 
r s 

b = 90 - ip t 



A 



l* t " x s !' if > x t " x s l - 180 

360 - |X t - X s | , if |A t - X s | > 180 



(2-7) 

(2-8) 

(2-9) 



1 . Point Azimuth Solution 

The value of arc a connecting the DF site and test 
point is given by the law of cosines for sides, 

cos a = cos b cos c + sin b sin c cos A (2-10) 
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There is no quadrant ambiguity since the cosine takes nega- 
tive values in the second quadrant and arc a is always less 
than or equal to 180 degrees because of the restriction on 
angle A to one hemisphere. Since the cosine of angle A is 
of interest in determining arc a, both the absolute value 
sign and the condition can be removed by applying a reduction 
formula. From equation (2-9) , 

cos | A — A 
t s 



cos A = 



cos (360 - \X - A g ), if |X fc - A | > 180 



But, through reduction, 

cos (360 - | X^_ - | ) = cos | A^_ - A | 



The absolute value sign can be dropped since cos (A^_ - Ag) = 
cos (Ag - A^_) . From this result and equation (2-10) 

cos a = cos b cos c + sin b sin c cos (A. - A ) (2-11) 

t s 

The vertex angle B at the DF site can be found by applying 
one of the laws of cosines for sides , 



cos B = 



cos b - cos c cos a 
sin c sin a 



( 2 - 12 ) 



Equation (2-12) is valid for test points either east or west 
of the DF site at vertex B. 

In order to make comparisons with a signal azimuth, 
angle B must be converted to an azimuth from true north, 
denoted 8. Visually, the conversion is based on whether the 
test point is eastward or westward of the DF site with res- 
pect to the proper pblar triangle. Mathematically, the 
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conversion is based on the magnitude of the longitudes of 
the DF site and the test point, but must also account for 
the discontinuity at the prime meridian. There are four 
general cases with two distinct outcomes as illustrated in 
Figure 4 . 

2 . Point Azimuth and Signal Azimuth Difference 

From a given DF site, the difference between the 
azimuth of any test point and a signal azimuth is the smaller 
of the two possible outcomes, depending upon whether the 
difference is measured clockwise or counter-clockwise. The 
azimuth difference is denoted w, and the azimuth discontinuity 
at true north requires that the computation be conditional. 
With a clockwise rotation for the azimuths, the smaller 
angle of azimuth difference will be positive if 6 is 
leading and negative if 3 trails 0. The possible outcomes 
and their conditions are as follows. 

a. If | 6 - 6| < 180, then w = 8 - 0 

b. If | 6 - 0 | > 180 and 6 > 0, then w = 6 - 0 - 360 

c. If | 6 - 0| > 180 and 6 < 0, then w=|3-0 + 36O 

(2-14) 

The sets of conditions on computing 6 (B) in equation 
(2-13) and in computing w(6,0) in equation (2-14) yield six 
possible outcomes for computing azimuth difference as a 
function of the vertex angle B and signal azimuth 0. 
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Azimuth and Vertex Angle 



360 - B 



FIGURE 4 
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b - e 



B - 36 0 - 0 



0 ) = 



B + 360 - 0 
- B + 360 - 0 



- B - 0 

- B + 720 - 0 



Two utility variables may be introduced in order to express 
azimuth difference in a constant form with values dependent 
upon the existing conditions. 

w - pB + q - 0 (2-15) 



The pair of utility variables can take on any one of the 
following sets of values 



( 1 , 0 ) 



(1,-360) 



(p,q) = < 



(1,360) 

(-1,360) 



(- 1 , 0 ) 



(-1,720) 



Concluding development, angle B can be expressed as a differ- 
entiable function of sines and cosines of known values of 
longitude and geocentric latitude. From equations (2-11) 
and (2-12) and the application of an angle-difference rela- 
tion. 
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B = cos 



/ L 



sm c 



|cos b - cos c [cos b cos c + sin b si 

(cos X. cos X + sin X, sin X )"| 
t s t s J 

sin ^cos ■*" [ 

+ sin X sin X )] 



in c 



cos b cos c + sin b sin c (cos X. cos X 

t s 



Substituting the values for arcs b and c from equations (2-7) 
and (2-8) and using the reduction formulas cos (90 - x) = 
sin x and sin (90 - x) = cos x, 



B = cos 



^sin ip t - sin ip g [sin ip t sin ip s + cos ip t 
cos ip (cos X. cos X + sin X sin X )] 
sin ^cos [sin ip^ si 

t sln Vl) 



cos ip 



in ip + cos ip t cos ip s (cos X^ 



cos X + sin X 
s 



By multiplying the term sin ip through the interior brackets 

of the numerator, writing the first two resulting terms as 

2 2 
sin ip. (1 - sin ip ) , which is equal to sin ip cos ip , 

"E S t S 

cos ip can be cancelled from both numerator and denominator, 
s 



B = cos 



1 ^[sin 



cos ip - sin ip cos ip. (cos X. 
s T s t t 



cos X g + sin X^ sin X g ) 



| j sin ^cos 1 j^sin ip t 



sin ip + cos ip. cos ip (cos X. cos X 
s t s t s 



+ sin X. sin X ) 
t s J 



(2-16) 
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III. THE FIX SOLUTION 



In Figure 5, each of four DF sites have measured azimuths 
of a signal from the same emitter at an unknown location. 

The signal azimuths are plotted for illustration as dotted 
arcs. A test point, , is established, and azimuths 

from each DF site to P are drawn in solid arcs. A better 
fit can be found by moving P somewhere within the region- 
bounded by the dotted arcs . 

Given a number n of DF sites, each reporting the azimuth 
of a received signal with respect to its own location to a 
net control, the fix problem can be viewed as finding the 




Azimuth Differences 
FIGURE 5 
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point on the globe which minimizes the sum of some function 
of the azimuth differences co^ , for i = l,2,...,n. 

A. OBJECTIVE FUNCTION 

Since the azimuth differences are signed, minimizing 
their sum may produce a solution that lies outside the region 
bounded by arcs corresponding to signal azimuths when large 
positive and large negative azimuth differences cancel each 
other out. Minimizing the sum of the absolute values of 
azimuth differences avoids this problem of positive and 
negative differences; however, since the partial derivatives 
of absolute value are not computationally tractable, solution 
techniques are limited. Minimizing the sum of the squared 
azimuth differences provides an objective function that is 
differentiable and reflects absolute magnitude. Therefore, 
the solution of the fix problem is taken to be that pair of 
coordinates (X^, ip which will 

n 2 

Minimize E to. (X, ,ip.) (3-1) 

i=l 1 t t 

B. SECOND-ORDER GRADIENT SEARCH 

Dropping the t subscript' to simplify exposition, the 
objective function, equation (3-1) , can be represented as 
f (X,iJO . Given an initial test point (X ,\p) , a truncated 
Taylor's series for a function of two variables can be used 
to approximate the functional values of other test points 
(X + AX, \p + Ai/i) . 
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f ( A + AA , ip + Aip) = f(X,ip) + A\p + AA 



+ 1 (AM 2 a 2 f(X,») + AX „ aifU J I + 1 (W 2 Ai^li 

2 3 A 2 dAdlp 2 dip 2 



Expressing the approximation in matrix form and letting 
superscript t indicate the transpose. 



fX + A A' 



^ + Alp 



= f ( A , ^ + (AX, Aip) 



+ | (A A, A ip) 




Letting denote an initial test point (X,ip) , letting R 

denote the vector (AX,Aip)^, letting V denote the gradient 
- 2 

of f(Pj) and letting V denote the Hessian matrix of f at 
Pj, the approximation can be expressed as 

f (P . + R) i f(P.) + R t Vf(P.) + - R t V 2 f(P.) R 
3 3 3 2 3 

Then , 



f (P . + R) - f(P.) = R t Vf(P.) + ~ R t V 2 f(P.) R (3-2) 

3 3 3 2 3 



Starting from an initial test point at P ^ , the search pro- 
cedure can be understood as a technique to find a vector R 
such that f (P j + R) is the minimum of the objective function. 
Such a vector R also minimizes the left hand side of equation 
(3-2) . Since the quadratic form on the right hand side is an 
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approximation of the value of the left hand side, it must 
be near a minimum value when the left hand side is minimized. 
The approximation works both ways, so that a vector R which 
minimizes the right hand side must approximately minimize 
the left hand side. The right hand side is in a form that 
yields to analysis of its extreme values with respect to R 
at a fixed P ^ . One of the conditions necessary for the 
quadratic form to be at a minimum is that its gradient be 
equal to zero. Therefore, it is necessary to find a vector 
R which fulfills this condition. Taking the gradient of 
equation (3-2) with respect to R, 

R fc Vf (P.) + ^ R fc V 2 f (P.) R 
3 2 3 





— t — 




*1 -t 2 ,- , - 


V R 


R Vf P . 

L 3 J 


+ V R 


j R Z V y f (P j ) R 



= Vf(P.) + \ V 2 f(P.) R + - R t V 2 f(P.) 

J Z J Z J 

= Vf (P . ) + v 2 f (P . ) R 
3 3 

Setting the gradient of the quadratic form at R equal to zero, 

.. Vf (P .) + V 2 f (P .) R = 0 
3 3 

Then , 




V 2 f (P .) R = -Vf (P .) 
3 3 



2 

And, providing that V f(Pj) is invertible, 



R = - 



V 2 f (P 



.•] 



-1 



Vf (P.) 



(3-3) 
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R specifies both a direction and magnitude in which to 
move toward an extreme value of the objective function. 

This consequence of the Taylor approximation can be improved 
by applying a scalar coefficient, denoted <5, to R in order 
to adjust the magnitude and, if R points to a maximum value, 
the direction. The next test point is then given by 

P,,. = P . + 6R (3-4) 

3 + 1 3 

Combining equations (3-3) and (3-4) , 

P j+1 = Pj - 6[v 2 f(P.)] _1 Vf(P j ) (3-5) 



The search procedure can be terminated when further 
improvements are of little consequence, most readily deter- 
mined by comparing the distance between locations on succes- 
sive iterations. 

In order to solve equation (3-5) , it is necessary to 
compute the inverse Hessian of the objective function. 



V 2 f (P) 



1-1 




3 2 f (P) 
3X3t|> 

a 2 f ( p) 

3X 2 



where the determinant of the Hessian 



det [v 2 f (P)l = — -j-J- 


3 2 f (P) 


3 2 f (P) 


1 J sx 2 


3^ 2 


dXdip 



Substituting the inverse Hessian in equation (3-5) and expand- 
ing the equation by matrix multiplication of the inverse 
Hessian times the gradient yields a new test point with 6 the 
only unknown. 
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p 



j+1 




The solution of equation (3-6) in turn requires the cal- 
culation of all first and second partial derivatives of the 
objective function. 



3f (P) 

d\ 



n 2 1 
E w 7 

i=l 1 



3A 



n 

E 2p . (p . B . + q . - 0 . ) 
i=l 11 1 1 i 



3B. 

l 

3 A 



(3-7) 



Similarly , 



3f (P) = 

dip 



n 



3B . 

l 



E 2p.(p.B. + q. - 0.) — 
i i i l l dip 



(3-8) 



3 2 f (P) 
3 A 2 



3 2 f (P) 



dlp‘ 



S 2 f (P) 

dXdlp 



n 

= E 
i=l 



n 

= E 
i=l 



n 

= E 
i=l 



3 B . 



3 2 B. n 



2p i (stH + 2p i (p i B i + 91- «i) 



(3-9) 



o i 9B. 

2p i 5T 



\ 2 



3 2 B . 

+ 2p . (p.B. + q. - 0.) — 

r l *1 l l . ,2 

3 ip 



(3-10) 
3 2 B. 



_ 3 B . 3 B . _ _ 

2PT 5T 1 57TT-+ 2 P^ (P,- B ; + q,- " 9,*) 



i 3 A 3ip 



i *i i M i l 3 A3^ 

(3-11) 
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Solving equations (3-7) through (3-11) requires the 
lengthy procedure of taking all first and second partial 
derivatives of equation (2-16) , the expression for angle B. 
Solutions to these equations are contained in the FORTRAN 
Fix Program included in this paper. 



C. GOLDEN SECTION SEARCH 

From equations (2-15) and (3-1) , the objective function 
can be written, 

n 2 

Minimize E (p . B . + q. - 0.) , (3-12) 

. . 1 l l l ' 

1=1 



the sum of the squared azimuth differences. 

Letting x denote the indicated move in longitude and y 
denote the indicated move in latitude given by equation (3-6) , 



R = (AX, Aip) 
= (x,y) t 



Then angle B^ , given in equation (2-16), becomes at 



B^ = cos 



-1 



[• 



sin (^ j - <5y) cos^ - sin^ cos (tJr - <5y) 



cos (X . -<Sx) cosX . + sin (X . - 6x) sinX. 
j l 3 l 



• r _i 

in cos 

|sin (ip j-6y) sini|R + cos (i|r - <5y) cos^ jcos (X_. - Sx) 



cosX^ + sin(X_. - 6x) sinX 



(3-13) 



where subscript j denotes test point coordinates and sub- 
script i denotes DF site coordinates. Since test point 
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coordinates are arbitrarily known and DF site coordinates 
are constant for the problem, the substitution of equation 
(3-13) into equation (3-12) , the objective function, results 
in the problem of minimizing a function with a single unknown, 
6, the scalar coefficient. 



2 

Minimize f(6) = £ to- (6) 

i=l 1 



(3-14) 



The objective function could be evaluated at incremental 
values of 6, say over the interval (a,b) . Then that 6 which 
yields the smallest functional value would be chosen as the 
scalar coefficient to adjust the move indicated by the 
second-order gradient search. The Golden Section Search 
(reference 1) provides a fast and methodical approach to find 
this minimizing 6. Figure 6 illustrates the problem; however, 
the actual curve representing f(6) is unknown at the outset 
of the problem. 
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Based on the assumption that f ( 5 ) has a local minimum on 
(a,b) , the objective is to find that minimum without evalu- 
ating the function at too many points in (a,b). The Golden 
Section Search systematically eliminates portions of the 
interval by comparing functional values of pairs of points 
in (a,b) . From Figure 6, l ess than f( 6 ^); there- 

fore, the portion (a, 6 ^) can be eliminated in searching for 
the minimum, reducing the search interval to ( 6 ^,b). A new 
point then is selected for comparison with resulting 

in a second elimination of whichever end interval corresponds 
to the higher functional value. 

Letting I represent the length of the n ^* 1 search inter- 
val , 

I 1 = b - a 

Two different values of <5 are tested at each trial. By 
letting the two values be symmetric on I , the length of 
I n+ ^ is the same regardless of which end interval is dropped. 
Letting a denote the relative position of 6 ^ in 1^, 

0 = ( 6 1 - a ) / (b - a) 

Then 

6 1 = a (b - a) + a (3-15) 

Since 62 is symmetric with 6 ^, the relative position of 62 in 
I 1 is given by (1 - c) , so that 

1 - o = (6 j ” a )/ ( b “ a) / 
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and 



6 2 = (1 - o) (b - a) + a (3-16) 

Provided that 6 2 yielded the minimum functional value, as in 
Figure 6, then <$ 2 becomes the first element of the symmetric 
pair in I 2 and its relative position in I 2 is given by 

0 = (6 2 - 6 1 )/(1 - a) (b - a) 

= [(1 - c) (b - a) + a - 6j/(l - a) (b - a) 
from equation (3-16) , so that 

6 1 = (1 - a ) 2 (b - a) + a (3-17) 

From equations (3-15) and (3-17) , it follows that 

a (b - a) + a = (1 - o) 2 (b - a) + a 

a = (1 - a) 2 (3-18) 

If, on the other hand, 6^ yields the minimum functional value, 
then 6^ becomes the second element of the symmetric test pair 
in I 2 , and its relative position in I 2 is given by 

1 - o = (6 1 - a) / (1 - a) (b - a) 

So that 

6 1 = (1 - a) 2 (b - a) + a, 

as in equation (3-17) , and leading to the same result as in 
equation (3-18) . 

From equation (3-18) 

0 2 - 3a + 1 = 0 

which can be solved by the quadratic formula, so that 
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a = [3 ± (9 - 4)^/2 
= 0.382 

Substituting this result in equations (3-15) and (3-16), 

6 1 = 0.382 (b - a) + a (3-19) 

6 2 = 0.618 (b - a) + a (3-20) 

Generalizing, let b^ and a^ represent the endpoints of the 
th 

n 1 " search interval. Then 

6. = 0.382 (b - a ) + a (3-21) 

1 n n n 

6 0 = 0.618 (b - a ) + a (3-22) 

2 n n n 

t 

After the first comparison on 1^, either 6 ^ or 62 i- s the 
survivor from the preceding iteration , so that a comparison 
can be made after evaluating only the new test point symme- 
tric with the survivor. The Golden Section Search terminates 
when b n - a^ < e , where e may be tailored to the magnitude of 
the indicated move. 

The value of 6 which minimizes f ( 6 ) over (a^,b^) is 
entered in equation (3-6) , resulting in a new test point 

p j+1 = Uj, ^j) t + 6 (X,y) t (3-23) 

= P . + 6 R 
3 

D. FORTRAN FIX PROGRAM 

The equations presented in this paper were coded in 
FORTRAN. The resulting program has been successfully 
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executed on an IBM 360/67 computer system in the W. R. Church 
Computer Center at the Naval Postgraduate School. 

1 . Program Description 

The program contains an identification number, longi- 
tude, latitude and name for each DF site in an imaginary 
network. Inputs to the program are DF site identification 
number and signal azimuth. After DF site names and coordi- 
nates are correlated with input site identification numbers, 
necessary trigonometric functions of site coordinates are 
computed and stored in an array. Searching for the minimum 
sum of azimuth differences begins from (0,0) with a Golden 
Section Search on longitude using (360,0) as the indicated 
move. The result is a minimizing longitude on the equator 
from which a Golden Section Search is conducted on geocentric 
latitude with indicated move (0,180). From this minimizing 
geocentric latitude, another Golden Section Search is con- 
ducted on longitude with (180,0) as the indicated move. 

Such alternating searches on longitude and geocentric lati- 
tude continue until- diminishing returns are indicated by a 
new position being less than 500 nautical miles from the 
previous minimizing position. The initial interval for the 
Golden Section Search is set at (-1,1) , the endpoints for 
equations (3-19) and (3-20). Subsequent endpoints are 
modified by the search results , producing values for equa- 
tions (3-21) and (3-22) . 

The program then executes second-order gradient 
searches refined by Golden Section Searches on indicated 
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moves until the most recent move is less than some input 
termination value, say one-tenth of a nautical mile. Indi- 
cated moves, given by equation (3-6) , are unbounded and 
might occasionally go off the coordinate system of Figure 1. 
Scaling any such wild moves to the coordinate system by 
modulo arithmetic destroys both the direction and magnitude 
of the vector. The indicated moves are kept in bounds by 
conversion to a unit vector in radians until such moves are 
of magnitude less than 30 nautical miles, such a small step 
indicating that the solution is at hand. By converting the 
indicated move to a unit vector prior to refinement by a 
scalar coefficient from the Golden Section Search, moves 
from the second-order gradient search are limited to about 
3440 nautical miles per iteration. 

Search output on successive iterations is optional, 
depending on an operator entry. Unless search progress is 
requested, searching terminates after thirty iterations even 
if the latest step is greater than the input termination 
value. After every tenth second-order gradient search, if 
such a high number should be required, the search effort is 
jolted by conducting alternate Golden Section Searches on 
longitude with indicated move (90,0) and latitude with indi- 
cated move (0,90). Such Golden Section Search interruptions 
only yield a new test point when they reduce the current 
value of the objective function. 

Prior to solution output, geocentric latitude is 
converted to geographical latitude using equation (2-6), and 
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the solution is given in degrees and minutes of longitude, 
east or west, and latitude, north or south. 

2 . Results 

The program has worked successfully in over ten 
thousand runs against simulated targets spread randomly over 
the globe with various networks of DF sites and various 
numbers of sites intercepting a given target. In actual 
practice, prior information about a target's area of opera- 
tion might produce a better starting point than (0,0); how- 
ever, the time required to type such starting coordinates 
on a keyboard for input would exceed the total time required 
for the solution starting from (0,0), which averages less 
than a second and includes some diagnostic output useful in 
analyzing results. Figure 7 illustrates a representative 
set of search iterations with Golden Section Searches 
plotted as dashed lines and second-order gradient searches 
solid. 
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IV. TARGET LOCATION CONFIDENCE REGION 



That location yielding the minimum sum of squared azimuth 
differences, or least squares solution, would coincide with 
the location of the signal emitter if all measured signal 
azimuths with respect to the DF sites were free from error. 
Since then all great circle arcs corresponding to signal 
azimuths would intersect at a point, the minimum value of 
the objective function would be equal to zero, as is always 
the case when only two signal azimuths are available to 
compute a fix. Because the sites in a DF network do not 
produce errorless measurements , solutions rarely coincide 
with the emitter location, and there is the question of how 
close the emitter might be to the least squares solution. 

A. SIGNAL AZIMUTH ERROR DISTRIBUTIONS 

Supposing that signal azimuth measurements made at indi- 
vidual DF sites are not biased in any of the possible direc- 
tions, then it follows that signed errors have a mean value 
of zero. If such measurements were found to be biased in 
certain directions during routine calibration trials, then 
the amount of bias could be applied as a correction to 
subsequent measurements , producing an estimated mean value 
of zero. Lacking any reason to assume a different distribu- 
tion and following the lead given in the literature of 
direction finding, the error distribution can be taken to 
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2 

be normally distributed with mean zero and variance a/, 
assumed known and independent from site to site. 

The error distribution can be made to be symmetric 
around zero by following the sign convertion that clockwise 
errors are negative and counter-clockwise errors are posi- 
tive in sign. 



B. THE OBJECTIVE FUNCTION VALUE AS A CHI-SQUARE STATISTIC 



From standard textbooks on probability theory (for 

example, reference 2) , it can be verified that if a random 

variable is normally distributed with mean zero and variance 
2 

a , then the square of that random variable is governed by 

9 

the Gamma distribution with parameters (l/2,l/2a ). Further- 
more, the sum of n independent Gamma distributed random 
variables divided by the variance of the underlying normal 
distributions is a Chi-Square statistic with n degrees of 
freedom. Hence, where w denotes the azimuth difference, 
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(4-3) 



Since the normal distribution here referred to applies 
to the azimuth differences between the measured signal 
azimuth and the azimuth of the unknown emitter location, the 
uu refer to actual and unknown errors rather than calculated 
differences between solution azimuths and signal azimuths. 



36 



Equation (4-3) therefore specifies the distribution of 
values of the objective function (equation (3-12) ) at unknown 
emitter locations under the assumption that azimuth errors 
are normally distributed with known variance. 

C. CONFIDENCE LEVELS AND REGION BOUNDARIES 

Since the sum of squared azimuth errors divided by 
variances is a Chi-Square statistic under the previous 
assumptions, there is a specific probability, denoted a, 
that this sum will be less than or equal to any given value. 
Equivalently, for a fixed probability a, or confidence 
level, there corresponds an upper bound, denoted C ^ , 

which depends upon the degrees of freedom as well as the 
confidence level. 



The choice of confidence level, a, is arbitrary, The 
corresponding value for C , . can be found in Chi-Square 



For a given fix problem, the value of the sum of least 
squares, denoted f*, is the minimum possible value of the 
Chi-Square statistic for the observed signal azimuths. 
Furthermore, by evaluating the objective function at suc- 
cessive points outward in any direction from the least 
squares solution, it can be seen that the farther away the 
point evaluated up to the antipode, the higher the objective 
function value. The location corresponding to f* can 




(4-4) 



tables . 
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therefore be regarded as a point in a region which contains 
all the possible emitter locations for which the value of 
the objective function is less than or equal to an upper 
bound C a j , given the set of observed signal azimuths. 
Such a region is bounded by all locations for which the 
objective function has a value equal to C, , , and it can 



be stated that, for the observed signal azimuths, there is a 
probability a that the bounded region will contain the 
unknown emitter location. 

Daniels and Vajda (1951) challenge the utility of such 
an approach on the grounds that at given confidence levels 
the corresponding region may not exist with positive 
probability if n is greater than two. Their observation is 
a consequence of the fact that except for n equal to two 
all of the great circle arcs corresponding to signal azimuths 
seldom intersect at a single point thereby yielding f* 
equal to zero. 

Prior to the availability of the high-speed digital 
computer, the practical problems of producing such a region 
were certainly more striking than theoretical objections. 

Even so, the Daniels/Va jda criticism was more than casual, 
which may indicate a requirement for additional arguments 
to justify the use of such a confidence region. 

Suppose that f* is equal to k for some least squares 
solution, where 




(4-5) 
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and the probability p may be determined by interpolating 
with the value k in a table of Chi-Squared values. Daniels 
and Vajda are concerned that if the confidence level a had 
been initially specified at some value less than p, then 
there exists no region of locations having objective function 
values less than or equal to k , k being the minimum possible 
value on the globe for the given signal azimuth observations. 
There are at least two responses to such an occurrence and 
the accompanying empty region objection: (1) Increase the 

confidence level whenever the region does not exist at the 
prior confidence level so that C is greater than k, and 

(2) Routinely specify a confidence level of sufficient mag- 
nitude to make the occurrence of an empty region highly 
improbable . 

Either of the above alternatives insures the production 
of a non-empty confidence region. There remains the question 
of what is represented by a possible contention inferable 
from equation (4-5) , namely that positive probability p is 
attached to a single point, the least squares solution. Any 
such statement would be puzzling when approached with the 
viewpoint that specific random variables on a continuum have 
zero probability and that probability statements can only be 
made about boundaries for continuous random variables such 
as location coordinates. 

An answer may be attempted by first pointing out that 
azimuth errors are unknown and independent random variables, 
but they are attached to actual and unknown emitter azimuths 



39 



with respect to the set of intercepting DF sites. Observed 
signal azimuths have two components, the first of which, 
emitter azimuth, is dependent upon emitter location. The 
second component, azimuth error, is independent from site to 
site and is the theoretical basis for the Chi-Square contour, 
confidence region. Prior to the production of a particular 
confidence region, the normality assumption permits the 
specification of upper bounds for the Chi-Square statistic 
corresponding to various levels of confidence. During 
production of the particular confidence region, the require- 
ment that a single location must satisfy the observations 
places a lower bound on the Chi-Square statistic. 

Each emitter azimuth can depart from the observed signal 
azimuth, but there are conceptual bands (based .on the error 
distribution assumptions) bounding such departures for any 
given probability. When more than two such bands are over- 
lapped in the solution region, then a proportion of the 
possible combinations of such departures are ruled out by 
the constraint that azimuth errors are attached to emitter 
azimuths. Whenever f* is greater than zero, then the prob- 
ability, denoted p* , associated with f* as an upper boundary 
from equation (4-5) represents the proportion of impossible 
sums of squares out of those permitted by the density func- 
tions governing azimuth errors. 

A thought experiment may clarify these remarks. Suppose 
it were possible to specify all combinations of azimuth 
errors permitted by their respective density functions prior 
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to finding the least squares solution. There exists a loca- 
tion satisfying one of the postulated combinations of azimuth 
errors for a specific set of signal azimuths only if applica- 
tion of each of the postulated azimuth errors in the combina- 
tion to its corresponding signal azimuth yields a corrected 
set of emitter azimuths whose corresponding great circle arcs 
intersect at a point. Few such postulated combinations in 
the infinite set can be satisfied since the point of inter- 
section for any two corrected emitter azimuths uniquely 
determines the required azimuth corrections for the rest of 
the observations. Furthermore, there cannot exist a satis- 
fying location for any postulated combination of azimuth 
errors having a sum of squares less than f* since no such 
corrections are of sufficient magnitude to result in a 
common intersection. 

These suggestions can be illustrated by an extreme and 
improbable example of emptiness of the confidence region for 
three signal azimuth great circle arcs forming a cocked hat. 
Suppose that the relationships between the DF sites and the 
emitter are as given in Figure 8, and suppose that the errors 
associated with the three DF sites are all governed by the 
standard normal distribution (variance equal to one degree 
squared) . Suppose also that each of the three signal 
azimuths were in error by minus two degrees (where the 
probability that azimuth error is less than or equal to 
minus two degrees is equal to .0227) in order to form one 
of the two largest possible cocked hats for | | equal to 

two degrees. Then the least squares solution will be found 
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Bounded Signal Azimuths 



FIGURE 8 
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near the center of the cocked hat, and f* is about equal to 
twelve degrees squared. For three degrees of freedom, there 
was a .992 probability that the Chi-Square statistic would 
be less than or equal to twelve degrees squared, indicating 
an empty confidence region at the .992 level (or any lower 
levels). Now the probability that | | is greater than four 

degrees is equal to .000, suggesting that a boundary of four 
degrees on either side of each signal azimuth great circle 
arc should contain all possible emitter azimuths from which 
each signal azimuth represents a departure. The situation is 
illustrated in Figure 8. 

With each emitter azimuth free to be anywhere within its 
dashed band by the normality assumption, but dependent upon 
the actual emitter location and therefore constrained by the 
other dashed bands, the shaded intersection of the three 
bands must (for any non-negligible probability) contain the 
emitter. Moreover, the Chi-Square contour methodology would 
produce a circular confidence region within the shaded 
triangle for any confidence level greater than .992. The 
center of the shaded region corresponds to the least squares 
solution. Suppose a "center" to be the least squares solu- 
tion united with all points forming a neighborhood within 
radius e of the center. Then this "center" had a probability 
of .992 + f(e) of containing the emitter location. Thus 
there was a probability of .008 - f(e) of finding the emitter 
within a region bounded by a "centerless" circle inscribed 
in the shaded triangle. There is zero probability that the 
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emitter lies outside such an inscribed circle, or confidence 
region at the 100 percent level. The description of the 
confidence region as circular applies to Figure 8 only 
because the DF sites are equidistant and symmetric around 
the emitter, and because azimuth errors are minus two degrees. 

The implication of the foregoing is that if errors are 
normally distributed with known variance, then abnormally 
large cocked hats would be indicative of excellent solutions 
via the least squares approach. The point would be difficult 
for operators to accept since high confidence is normally 
associated with a set of great circle arcs corresponding to 
signal azimuths that nearly intersect at a point, and large 
departures are suspect. The discrepancy between the above 
conclusion and experience may be attributed to the fact that 
abnormally large cocked hats only occur when outliers (signal 
azimuth errors that do not obey the empirically estimated 
distribution) enter the data base. 

Sometimes strange relationships between the sizes of 
cocked hats and accuracy of locations (least squares solu- 
tions) can also be demonstrated with a network of three DF 
sites with errors governed by a Bernoulli distribution, 
having equal probability that azimuth error is plus or minus 
three degrees. Figure 9 illustrates the situation. Given 
an emitter at the center of the star, observed signal azi- 
muths are in error by plus or minus -three degrees, each 
with probability 1/2. All possible configurations of 
cocked hats are therefore represented on the figure for the 
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Bernoulli Random Error 
FIGURE 9 
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sketched relationship between DF sites and emitter location. 
Six of the possible cocked hats are the small triangles at 
the tips of the star. The other two possibilities are the 
union of the body of the star with the small triangles at 
points NE , SE and W, and the union of the body of the star 
with the small triangles at points E, SW and NW. The latter 
two cocked hats are abnormally large and have f* equal to 27 
degrees squared (which is over twice the value of a Chi-Square 
statistic at the .995 level for a normal error distribution). 
Yet, the center of these largest cocked hats coincides 
approximately with the emitter location, indicating excellent 
least squares solutions, and the centers of the triangles at 
the points of the star, the smaller cocked hats, are far 
from the emitter location, even though their f* value is 
relatively small. 

The probability of occurrence of such large cocked hats 
when error is governed by a normal rather than Bernoulli 
distribution is negligible, and they do not appear when data 
is simulated with a normal error distribution having reason- 
able variance. However, in about 400 simulations during 
development of the Fix Program they were forced to appear 
by governing error with such a Bernoulli distribution. As 
would be expected from Figure 9, a plot of location error 
versus f* verified that location error decreases as f* 
approaches its maximum of 27 degrees squared for this unusual 
error distribution. Such a simulation naturally yields about 
25 percent of the abnormally large cocked hats and 75 per- 
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cent triangles at the points of the star. With a random 
distribution of DF sites and signal emitters over the globe, 
the stars are not so well-shaped as that in Figure 9. 

If one regards the fix problem with least squares solu- 
tion and confidence region construction at confidence level 
a as a Bernoulli trial with success if the confidence region 
contains the emitter location and failure if the emitter is 
outside the confidence region, then the expected number of 
successes for m independent trials is am. That the Chi- 
Square contour confidence regions have normalized frequency 
of containment about equal to the probability of containment 
is therefore expected to be true and was demonstrated to be 
true by simulation. At confidence level a, on about a per- 
cent of the trials non-empty confidence regions are produced 
which contain the actual emitter. On the remaining Cl - a) 
percent of the trials , outcomes are distributed among empty 
regions and non-empty regions with the actual emitter out- 
side the confidence region. The lower the confidence level 
the higher the proportion of empty regions amongthe failures. 

D. FORTRAN CONFIDENCE REGION PROGRAM 

A plot subroutine is included in the Fix Program of this 
paper. Sample plots are included in the Computer Output 
section of the paper. The plot is a square grid of constant 
dimension, scaled in nautical miles to contain the confidence 
region. An asterisk on the grid indicates that the corres- 
ponding point has an objective function value greater than 
the Chi-Square boundary value for the specified confidence 
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level. The least squares solution is at the center of the 
grid and is indicated by the letter S. During program 
development, it was found that a 90 percent confidence level 
is sufficient to preclude the occurrence of empty confidence 
regions, and the sample regions are of that size. 

1. Program Description 

Given in a data statement of the subroutine are 
Chi-Square values corresponding to the numbers of possible 
degrees of freedom, up to n for a network of n DF sites, and 
the selected confidence level. The confidence region boun- 
dary value is selected according to the number of signal 
azimuths available to find a solution. 

The dimension of the grid is 41 rows and 51 columns, 
tailored to the line printer distances between rows and 
columns, so that the ratio of the distances between rows and 
columns is 4:5. Depending upon the scale of the grid, coor- 
dinates on the globe are calculated to match each of the 
2091 points on the grid. Curvature of the earth calculations 
are used to determine distances between points on the grid, 
thereby avoiding distortion inherent in planar assumptions. 

At each of the points on a row, the sum of squared 
azimuth differences is computed for comparison with C a )/‘ 
an asterisk being assigned for grid entry if C ^ is 
exceeded. After the confidence region is exited on a given 
row, remaining points on the row are assigned asterisks 
without computation. Similarly, after the confidence region 
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is completed, the rest of the grid is filled with asterisks 
without computation. 

The area of the confidence region is estimated by 
counting the blanks and multiplying their sum by the area 
represented by each cell, 80 square miles in the case of 
a 400 by 400 nautical mile grid. If computations on the 
first row yield any points with objective function value less 
than or equal to C ^ ^ , then the grid is rescaled to cover 
a greater area, and the computation starts over. If the 
confidence region leaks out on either of the sides of the 
grid during computation, then again the grid is rescaled, 
and computation starts over. In practice, the initial area 
covered by the grid can be tailored to the network and 
expected size of confidence regions. 

2 . Results 

If the solution search in the Fix Program should 
terminate in the region of the least squares solution before 
the least squares point is found, then the confidence region 
is noticeably asymmetric about the center at S. The lack of 
symmetry indicates that the fix solution should be recomputed 
with a smaller termination step. The center of the confi- 
dence region, or least squares location, could also be 
estimated from the grid itself, but only to an accuracy 
limited by the dimensions of a cell. 

The confidence region subroutine was used successfully 
in over a thousand runs with simulated data. Its chief 
deficiency is the computer time required to produce the grid 
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entries, averaging slightly over ten seconds for numbers of 
signal azimuths per target varying uniformly from three to 
nine. On a slower machine than the IBM 360/67, the subrou- 
tine would require revision to be usable for production. 
Sampling every other or every third cell would be possible 
avenues for reducing computation time, but the present 
version was preferred for a more accurate estimate of con- 
fidence region area. A more sophisticated method for reduc- 
ing computation time would be to use a modified Golden 
Section Search to find the confidence region boundary outward 
from the center of the region, store coordinate pairs of the 
boundary in an m by n array, and then print the boundary on 
a scaled grid. 
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V. 



INTERCEPT GEOMETRY 



The Fix Program produces several statistics that are 
useful in describing intercept geometry. The first group 
of these statistics is related to signal azimuths and is 
used to plot position lines as linear approximations of the 
great circle arcs in the region of the solution. In the 
second group are some values to which location error is 
sensitive. 

A. PLOTTING POSITION LINES IN SOLUTION REGION 

Following the optional summary of search progress, the 
least squares solution and the value of f * , the computer 
output has a table of values in which each row corresponds 
to one of the intercepting DF sites. Row entries include 
the site identification number, the site name, the azimuth 
of the solution, the signal azimuth, their difference, the 
range in nautical miles from the solution to the site, and 
three values used to plot the position line relative to the 
least squares solution. These values are calculated by 
analogy with a right spherical triangle illustrated in 
Figure 10 . 

Figure 10 is labelled to correspond to the base of the 
polar triangle of Figure 3 with the DF site at vertex B and 
the least squares solution at Vertex C. Vertex angle A is 
a right angle formed by connecting C to the closest point of 
approach (CPA) of the great circle arc related to the signal 
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Right Spherical Triangle 
FIGURE 10 



azimuth 0. Arc a, connecting the DF site and least squares 
solution, is related to the solution azimuth, denoted 3*, 
and is given by equation (2-11) , using the polar triangle of 
Figure 3. By the lav; of sines for spherical triangles, 

sin b = sinto sin a (5-1) 

The curvature of earth distance covered by arc b is 
closely approximated by multiplying arc b in radians by the 
mean earth radius of 3439.9 nautical miles, yielding the 
nautical mile departure of the CPA from the solution. 

The site azimuth with respect to the solution, denoted 
y*, is found after 1 calculating angle C of Figure 3 in the 
same manner that vertex angle B and azimuth 3 were calculated 
in equations (2-12) and (2-13) . Then the azimuth of the CPA 
with respect to the solution is the sum of y* and c, in the 
case illustrated. Finally, the azimuth of the DF site with 
respect to the CPA is found by subtracting 90 degrees from 
the sum of y*and C, in the case illustrated. Similar argu- 
ments apply to three other general cases, solution east or 
west of the DF site and to positive or negative, and resulting 
conditional equations are included in the Fix Program. 






The azimuth and distance of the CPA with respect to the 
solution are the components of a vector which can be plotted 
on polar coordinate paper. A line of position perpendicular 
to this vector approximates the great circle arc in the 
region of the solution. The direction of the intercepting 
DF site is given by the azimuth value in the last column of 
the output row. The plot can also be constructed with a 
parallel ruler by orienting the ruler along the site azimuth 
with respect to the CPA, moving the free edge the distance 
to scale of the CPA toward the azimuth of the CPA. Such 
lines of position are plotted to the grid scale on several 
of the confidence region plots provided in the computer out- 
put section of this paper. 

B. MEAN ACUTE ANGLE OF INTERSECTION 

There are factors that may be said to reduce location 
error, given that signal azimuth error is present but un- 
known. One such factor is the number of DF sites providing 
signal azimuths to make a fix: the more independent measure- 

ments, the better the emitter location estimate by least 
squares solution. Another factor is the layout of the DF 
network, the obvious supposition being that a network spread 
over a large area produces better locations of globally 
distributed targets than a network confined to a small area. 

A slightly less obvious statistic is the magnitude of 
the acute angle of intersection of great circle arcs corres- 
ponding to signal azimuths. For two DF sites participating, 
the dependence of location error on the acute angle of 
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intersection is illustrated in Figure 11. In case a, two 
lines of position intersect at about 90 degrees; in case b, 
the acute angle of intersection is about ten degrees. The 
difference in the sizes of regions contained by the dashed 
lines demonstrates the effect of two degree uncertainty of 
measurement on either side of a signal azimuth on the two 
different intersections. 

Although the confidence regions described in Section IV 
of this paper do convey this dependence, that for case a 
being a near circle contained in the dashed lines and that 
for case b being an elongated ellipse of orientation and 
dimension similar to to the dashed region, the acute angle 
of intersection statistic was found to be significant enough 
to warrant isolated analysis. It can be demonstrated 
graphically that the dependence also exists for n greater 
than two. For all cases, the dependence can be examined by 
calculating the mean value of all the acute angles of inter- 
section formed by all possible permutations of two lines of 
position, given by (n^* The mean acute angle of inter- 
section is part of the Fix Program output. Dependence of 
location error on this statistic is examined in later sec- 
tions of this paper. For individual pairs, the acute angle 
of intersection is closely approximated by taking differences 
from the site azimuths with respect to the CPA. 
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b. 10 Degree Acute Angle of Intersection 



FIGURE 11 
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C. MEAN ARC STATISTIC 

A more subtle factor affecting location error is a 
function of arc length or distance from the DF site to the 
emitter location. The arc from the DF site to the emitter 
and all other arcs representing possible departures due to 
erroneous signal azimuth measurements converge at the DF 
site and at its antipode. Therefore, the maximum nautical 
mile departure of a line of position from the emitter loca- 
tion for constant azimuth error occurs at 90 degrees of arc, 
or about 5403 nautical miles, from the DF site. 

The arc from the DF site to the least squares solution 
is given by equation (2-11) . The arc statistic is equal to 
this arc, or to its supplement if the arc is greater than 
90 degrees. The dependence of location error on the arc 
statistic is illustrated in Figure 12, where the solid line 
represents an arc from the DF site through the emitter loca- 
tion to the DF site antipode, and the dashed lines represent 
three degree departures from this arc. 

The distance from the site-to-emitter arc to the depar- 
ture arcs increases non-linearly from zero at the DF site 




Site 



Antipode 



Relating Error to Emitter Distance' 
FIGURE 12 
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and its antipode to 180 nautical miles at the midway point. 
In a later section of this paper, location error sensitivity 
to the mean arc statistic for all DF sites participating in 
the location is examined. 
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VI. SIMULATION PROCEDURE 



A modified version of the Fix Program was used to simu- 
late operation against targets. Targets were generated by 
a subroutine which provided longitude and geocentric lati- 
tude describing random locations uniformly distributed over 
the globe. Such coordinates were saved for later comparison 
with the least squares solution. 

The number of intercepting sites for each target was 
either specified for simulations on fixed degrees of freedom, 
or randomly generated according to a uniform distribution 
from three to nine. In both cases, the specific intercept- 
ing DF sites were randomly selected from a set of 19 DF sites 
composing the particular network being simulated. Azimuths 
of the target location with respect to the participating DF 
sites were then calculated. For each such emitter azimuth, 
a random error was added, with direction and magnitude of 
the error governed by the standard normal distribution. The 
emitter azimuth with error was then treated as a measured 
signal azimuth for input with the site identification number 
to the Fix Program. 

During program de\elopment, the simulation procedure was 
used with zero signal azimuth error generated to insure that 
the least squares solution of the search procedure coincided 
with the target coordinates to at least an accuracy of three 
decimal places. A search termination value of one-tenth of 
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a nautical mile was found to be adequate, and the present 
combination of Golden Section Searches and second-order 
gradient searches was found to be satisfactory during the 
development simulations. The Fix Program has also been 
successfully tested with up to 19 participating DF sites on 
a single target, with a set of DF sites contained in a 180 
by 240 nautical mile area, and with higher variance error 
distributions, including outliers. 

A. GEOCENTRIC AND GEOGRAPHICAL LATITUDE 

A comparison of least squares solutions using geocentric 
and geographical latitude in Fix Program calculations was 
made in order to examine the accuracy gain of the geocentric 
model. Errorless signal azimuths of random targets were 
computed using geocentric latitude in calculations. A least 
squares solution was then found by using geographical lati- 
tude in calculations, and the conversion of the solution to 
geographical latitude was omitted. This distance between 
the geographical coordinates of the target and the pseudo- 
geographical coordinates of the least squares solution 
measures the departure from the least squares location due 
to use of geographical latitude on an assumed sphere. 

Radial departures for two degrees of freedom were occa- 
sionally devastating (424 nautical miles in one case). Using 
an imaginary network, called the national network, and degrees 
of freedom uniformly distributed from three to nine, the mean 
radial departure for 148 samples was 6.8 nautical miles. 

The direction of such departures tends toward the equator for 
targets in both northern and southern hemispheres. 
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VII. LOCATION ERROR ANALYSIS WITH FIXED DEGREES OF FREEDOM 



The sensitivity of location error to several factors was 
examined in simulations with degrees of freedom fixed from 
two to seven. From 100 to 200 samples were generated at each 
level for each of the networks. All of the following plots 
are drawn to the same scale to simplify between plot compari- 
sons of the effects of different degrees of freedom and 
different network configurations. Values too large for the 
common scale are entered in a row at the upper limit of the 
plot above an off of scale indicator. 

A. COMPARISON OF ERROR ON NATIONAL AND GLOBAL NETWORKS 

Member sites of the national network varied in longitude 
from Attu to Boston and in latitude from Anchorage to Hono- 
lulu. The specific sites are identified in the second data 
statement of the Fix Program. Sites in the widely distri- 
buted global network included Anchorage, Bombay, Boston, 

Buenos Aires, Cape Town, Christchurch, Dakar, Honolulu, Houston, 
Lima, Naples, Papeete, Perth, Recife, Seattle, Seychelles, 
Tehran, Tokyo and Wake. 

Figure 13 is a plot of mean location error from each of 
the sample groups against the number of intercepting sites. 
Dashed lines connect the means of each network. The plot 
shows nearly double the location error in the national net- 
work even at higher degrees of freedom. 
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With all other factors random, the average accuracy of 
a network is dependent upon the network layout and the 
variance of signal azimuth measurement error. For a given 
network, the mean location error would vary in the same 
direction as changes in such variance, in all of the follow- 
ing simulations held constant at one degree. 

B. ERROR VERSUS MEAN ACUTE ANGLE OF INTERSECTION 

Based on samples from the national network, Figures 14, 

15 and 16 are scatter plots of location error versus the 
mean acute angle of intersection explained in section V-B. 
Degrees of freedom were fixed at three, five and seven, 
respectively . 

The strong correlation in the national network can be 
interpreted as a reduction of both the mean and variance of 
location error as intersections of great circle arcs are less 
acute. Rough scatter plots for two, four and six degrees of 
freedom showed the same correlation, naturally most pro- 
nounced at two degrees of freedom. The set of plots shows 
a shrinking of the range of the mean acute angle of inter- 
section statistics as degrees of freedom are increased. 

Figures 17, 18 and 19 are similar scatter plots with 
entries from simulations on the global network. The corre- 
lation is still evident at three degrees of freedom but be- 
comes masked at higher levels as the layout of the network 
reduces the chance of producing a combination of target and 
intercepting sites with a small mean acute angle of inter- 
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section. Consequently, shrinking of the range of the 
statistic is extremely evident on Figures 18 and 19. 

C. ERROR VERSUS MEAN ARC STATISTIC 

Figures 20, 21 and 22 are scatter plots of national 
network location error versus the mean arc statistic 
(explained in section V.C) for three, five and seven degrees 
of freedom, respectively. At all levels, the large location 
errors tend to be associated with the high end of the range 
of mean arc statistics. The strength of the correlation 
appears to increase as de'grees of freedom increase on 
successive plots. 

Figures 23, 24 and 25 are similar plots for the global 
network. Again, the correlation is masked after three 
degrees of freedom with the shrinkage effect extremely pro- 
nounced as the degrees of freedom increase. 
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IIX. LOCATION ERROR ANALYSIS WITH VARYING DEGREES OF FREEDOM 



Simulations were conducted on both networks with degrees 
of freedom generated uniformly from three to nine. Chi- 
square contour confidence regions were generated for each 
target's least squares solution. Using a confidence level 
of 90 percent, national network confidence regions contained 
90.38 percent of the generated targets. Results for the 
global network were 90.44 percent. Such close agreement 
between the confidence level and containment in the simula- 
tions is dependent upon a reliable standard normal azimuth 
error generator. In production use, such reliability could 
only be achieved by having a good estimate of DF site 
variances. 

A. NATIONAL NETWORK 

Figure 26 is a scatter plot of confidence region area 
versus location error. Again using a common scale to aid 
network comparison, eleven points off the scale are entered 
above or to the right of an indicator. Although there is 
one hopeful cluster near the plot origin, the large variance 
in location errors at given confidence region areas is 
indicative of the futility of predicting location error 
based on the size of this reliable confidence region model. 

Figure 27 is a multi-level histogram of location error 
at the various degrees of freedom, demonstrating the same 
sort of correlation evident in Figure 13. Figure 28 
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illustrates the correlation between location error and the 
mean acute angle of intersection, present even with mixed 
degrees of freedom. The same holds true in Figure 29, a 
plot of location error versus the mean arc statistic. 

B. GLOBAL NETWORK 

Figures 30 through 33 are similar plots for the global 
network. On Figure 30, location error and confidence region 
area appear independent. Figure 31 shows good correlation 
between location error and degrees of freedom. Figures 32 
and 33 show nearly complete masking of the previously 
established correlations due to the layout of the network. 
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IX. CONCLUSION 



The Fix Program and Chi-Square contour confidence region 
subroutine could be integrated into the software of position 
location systems by direction finding. The present program 
(with target and azimuth error generators included) could be 
used to simulate an existing DF network for accuracy studies 
on potential targets in different regions of the globe or for 
studies on optimal placement of DF sites in a proposed net- 
work . 

The relationships demonstrated in the sections on loca- 
tion error analysis indicate that larger location errors 
accompany contraction of the network layout; however, this 
factor would be mitigated by a practice of assigning regional 
rather than globally distributed targets to the network for 
prosecution. 

Countering the widespread belief that small confidence 
regions indicate accurate locations and large confidence 
regions indicate poor locations, it has been demonstrated 
that for a reliable model location error and confidence 
region area are nearly independent. 

The relationships examined as parts of location error 
analysis suggest an alternative answer for the question of 
how accurate is a given location. For an existing network, 
data could be collected from locations of known transmitters. 
Such data would include statistics on number of signal 
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azimuths available to make the fix, mean acute angle of inter- 
section and mean arc statistic. With known location error 
as the dependent variable, multiple linear regression 
analysis could be used to develop regression and correlation 
coefficients for predicting location error as a function of 
independent variables measured or calculated at the time of 
location . 

The workability of such an approach was not demonstrated 
in the present paper for several reasons: (1) The effort 

required to perform the multiple linear regression (using a 
standard program from the IBM scientific subroutine package) 
would be a project in itself and would be best reserved for 
actual data from an existing network. (2) System operators 
would be aware of other measurable factors for inclusion as 
independent variables in the regression analysis, such as 
variance on given signal azimuth measurements, signal 
strength and measurement duration. 
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FIX PROGRAM 



RFAL LAMBDA,N,NL,NP,NLL,NPP,NPL 

DIMENSION DFSITE(6,19) ,D0MEGA(19) , SYMBOL ( 4 ), BETA! 19), 
1STAA7 ( 19 ) 

COMMON/ SITES/S I TE ( 10 , 19 ) , NUMBER 
COMMON /SIGNAL/ THETA (19) 

DATA T WOP I /6. 283 18 5 /,P 1/3. 14 1593/, RADDEG/ 57. 29578/, 

I DEGRAD/1 .745329E—2/ , FNC ECC/ 1 .00 6768/ , SCREEN/ .9999999/, 



2PIHALF/1 .57079 


6/ , SYMB CL/ ' N ' , ' 


S 


' , ' E« , 


• W' / 




DATA 


DFS ITE/ 














IDENT 


LONG 


LAT 


SITE NAME 




1 


1 ., 


210. 117, 


61.217, 


* ANCH' 


t 


' ORAG' 


, ' E 


2 


2. , 


173.250, 


52.933, 


' ATTU ' 


* 


i t 


, ' 


3 


3. , 


288. 933, 


42.350, 


' BOST' 


* 


' ON * 


, ' 


4 


4., 


262.600, 


27. 800, 


•CORP • 


t 


'US C' 


, 'HRI S 


5 


5. , 


254.983, 


39.717 , 


« DENV ' 


f 


« ER • 


, ' 


6 


6. , 


276.950, 


42.333 , 


• DETR« 


f 


' 01 T ' 




7 


7., 


253.517, 


31.750, 


•EL P ' 


t 


'ASO ' 


, ' 


8 


8. , 


247.983, 


46.6 00 , 


• HELE • 


* 


' NA « 


, ' 


9 


9. . 


202. 133, 


21.317, 


' MONO' 


t 


' LULU' 




1 


10 . , 


241.750, 


34.050, 


•LOS ' 


1 


• ANGE • 


, ' LES 


2 


11. , 


279. 800, 


25.767 , 


• Ml AM* 


f 


' I • 


, ' 


3 


12., 


182. 633, 


28. 217, 


•MI DW' 


1 


5 AY ' 


, ' 


4 


13. , 


266.783, 


44.983 , 


• M I NN ' 


* 


« EAPO • 


, 'LIS 


5 


14. , 


273. 200, 


36.150, 


' NASH' 


t 


' VILL ‘ 


, ' E 


6 


15., 


269.383, 


29. 967, 


'NEW • 


t 


' ORLE • 


, ' ANS 


7 


16. , 


283.767, 


38 . 667 , 


' NORF « 


f 


'OLK ' 


, ' 


8 


17. , 


240. 200, 


39, 517 , 


• RENO* 




, t 


, ' 


9 


18 


237.667, 


47.600, 


'SEAT' 


f 


‘ TLE ' 


, ' 


1 


19. , 


262.667, 


37.683, 


' W1CH' 


9 


' ITA • 


, ' 



0900 



1000 



, 

/ 

READ (5,1) KILL, CUT 
READ (5,1) NUMBER 
IF ( NUMBER. EO. 0) GO TO 4200 
WRITE (6,9) 

FORMAT PI') 

FORMAT (I2.F8.3) 

DO 1000 K = 1, NUMBER 
R F AD (5,1 ) ID, THETA (K) 

THETA ( K ) = THETA ( K)*DEGRAD 
SITE(1,K)=DFSITE(2, ID)*DEGRAD 
S I TE ( 2 , K) = SIN( S ITE( 1 ,K ) ) 

SITF(3,K)=C0S( S I TE ( 1,K) ) 

SITE(4,K )=ATAN ( 0 .99 32773+TAN ( D FS I TE ( 3,1 D)*DEGRAD) ) 
SITE(5.K)=SIN(SITE(4,K) ) 

SITE! 6 , K ) =C0 S( SITE( A , K ) ) 

SITE(7,K)=DFSITE(1, ID) 

SITF(8.K)=DFSITE(4, ID) 

SITE(9,K ) = CF SI TE ( 5, ID) 

SITE(10,K)=DFS I T E ( 6 , I D ) 

C PERFORM ALTERNATING GOLDEN SECTION SEARCHES CN LONGITUDE 
C AND LATITUDE 
PSI=0. 

LAMBDA = 0 . 

S AV EL = L AMBDA 
S AVFP = P S I 
DEGF = 6 1 5600. 

CO 1400 J = 1, 5 
DI V = J 

X=l. /DT V* TWO PI 
Y = 0. 

CALL SEARCH! LAMBDA, PS I , X, Y, DEGF ) 

IF (KILL. EG. 2) GO TC 1100 

CALL SH I FT ( SAVEL , SA VEP , LAMBD A , P S I , 1 ,DEGF ,D I ST) 
CONTINUE 
X = 0. 

Y=1./DIV*PI 

CALL SEARCH! LAM BOA, PS I,X,Y,DEGF) 

IF (KILL. EG. 2) GO TC 1200 

CALL SHIFT ( SAV EL , SAVEP , LAMBDA ,PSI , 1 , DEGF , DI ST) 

GO TO 1300 

CALL AZDIST( SAVEL, SAVEP, LAMBDA, PS I ,AZGO, CIST) 



1100 



1200 
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SAVEL=L AMBDA 
SAVEP=PSI 

1300 IF (DI ST.LE.500. ) GO TO 1500 
1400 CONTINUE 
1500 P=0 

1600 M = M + 1 

COSLAM=CCS( LAMBDA) 

SINLAM=SIN( LAMBDA) 

COSPSI^COSIPSI ) 

S I NPS I = S I N ( PS I ) 

FL = 0 . 

FP=0 , 

FLL -0 . 

FPP=0, 

FPL=0 . 

C EVALUATE FIRST AND SECOND PART I ALS OF OBJECTIVE FUNCTION 
DO 2450 K = 1, NUMBER 
A= 1 . 

C=0. 

UA = S I N P S 1 ^ S I TE (6,K) 

U B = COS P S I*S TTE ( 5,K ) 

UC=COSLAM* SI TE { 3 , K )+S I NL AM*S IT E ( 2 , K ) 

UD=-SINLAM«SITE{ 3, K )+C0 SL AM* SI TE < 2 , K ) 
UE=C0SPSI*SITE(6,K) 

UF=SINPSI*SITE ( 5 * K ) 

UG=UF+UE*UC 
UH=U B— U A'->UC 
UI =-UE*UD 
N=UA-UB-nJC 
NL =— U 3*UD 
NP = UE-f-UF*UC 
NLL=UB*’JC 
NPP=-N 
NPL=UF*UD 

IF (UG.GE.l.) UG=SCREEN 
IF (UG.LE.-l.) UG=-SCREEN 
D=S IN ( A RCO S ( UG ) ) 

UJ=1 ,/D 
DL=UJ* UG*UI 
DP =— U J*-UG*UH 

CLL = UJ* (UG-*UE*UC+U I*UI-DL*DL ) 

DPP=U J -4 ( UG* UG- UH*UH-OP* DP ) 

CPL=UJ*( UG*UA*UD+UH*UI-DP*DL ) 

COS B=N*U J 

IF (COSB.GF.l.) COSB=SCREEN 
IE (CDSB.LE.-l. ) CO SB=- SCREEN 
E=ARCOS ( COSB ) 

COTB=COT AN ( B ) 

CSCB=1 ./SIN(B) 

BL=UJ* ( DL*COTB— NL*CSCB ) 

BP = UJ*< DP*COTB-NP*CSCB) 

UK =U J*CS C 3 

B L L = UK* ( LJ*OL* { NL-DL* COS B ) +DLL*COS 3-BL* DL*CS CB-NL L + 
1BL*NL*C0TB ) 

BPP=UK* (UJ*DP* (NP-DP*COSB)+CPP*COSB-BP*DP*CSCB-NPP+ 
1BP*NP*C0TB) 

BPL=UK*( UJ*DP*(N L-DL* COSB) +DPL*COSB-BP*DL* CSC B-NPL+ 

1 BP*NL*COTB) 

TE S T = S I TF ( 1 ,K) — LAMBDA 
IF (TEST) 1700,2000,1800 
1700 IF (TEST+PI) 1900,2000,2000 
1800 IF (TEST-PI) 1900, 2CC0, 2000 
1900 A=— 1 . 

C=TWOPI 

2000 TEST=A*B+C— THETA(K) 

BETACK ) = (A-*B+C JtRADCEG 
IF ( ABS ( TEST)- PI ) 2400,2400,2100 
2100 IF (TEST) 2200, 2400,2300 
2200 C = C+T WOP I 
GO TO 2400 
2300 C=C— T WOP I 
2400 CPEGA = A* E + C-THET A ( K ) 
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DC MEGA ( K) =OMEGA*RADDEG 
UL=2.*A*(A*B+C-THETA(K ) ) 

FL=FL+UL*BL 
FP = FP + UL ; rBP 

FLL = FL L +2 . *BL* 3L +UL "BL L 
FPP=FPP+2.*BP* BP + UL^BPP 
FPL=FPL+2.*BP*BL+UL*BPL 
2450 CONTINUE 

UM = 1. /( FI L rt FPP-«=PL*FPL ) 

X= ( F P L * F P— F P P s ' : F L )*UM 
Y=(FPL*FL-FLL*FP )*UM 

C IF EITHER X OR Y ARE GREATER THAN PI AND UNTIL WITHIN 30 
C NM OF THE INDICATED MOVE, CONVERT (X,Y) TO A RADIAN UNIT 
C VECTOR 

IF (AMAX1C X, Y) .GE. PI) GO TO 2500 

GOX=LAMBCA+X 

GCY=PSI + Y 

CALI AZDIST! LAMBDA , PS I , GOX ,GOY , AZGO , DI ST ) 

IF (DIST.LE.30. ) GO TO 2600 
2500 CONTINUE 

DIST=SORT(X*X+Y*Y) 

X = X/ D I ST 
Y=Y/DI ST 
2600 CONTINUE 

C REFINE INDICATED MOVE BY A GOLDEN SECTION SEARCH 
CALL SEARCH! LAMBDA, PS I , X , Y , D EG F ) 

IF (KILL.EQ.2) GO TO 2700 

CALL SHIFT! SAV EL , SAVEP, LAMBDA, PS 1 , 2, DEGF, DIST ) 

GO TO 2800 

2700 CALL AZDIST (SAVEL, SAVEP, LAMBDA, PSI, AZGO, DIST) 

SA VE L=L AMBDA 
SA VEP = P S I 

C CHECK FOR MOVE LESS THAN MINIMUM STEP 
2800 IF (DIST. LE. CUT) GO TO 3300 

C SEARCHING STOPS AFTER 30 ITERATIONS UNLESS SEARCH PROGRESS 
C IS REQUESTED 

IF (KILL. EC. 2) GC TO 2900 
IF (M.EQ.30) GO TO 3300 
2900 CONTINUE 

C IF THF SOLUTION HAS NOT BEEN FOUND, APPLY THE GCLDEN 
C SECTION SEARCH AFTER EVERY TENTH ITERATION. 

IF ( NOD ( M, 10) . NE.O ) GO TO 1600 
X = P I HAL F 

Y = 0 . 

CALL SEARCH! LAMBDA, PSI ,X,Y,DEGF) 

IF (KILL. EC. 2) GO TO 3000 

CALL SHI FT (S AV EL, SAVEP, LAMBDA, PSI, 1,DEGF, DIST) 

GO TO 3100 
3000 SAVEL=L AMBDA 
SAVEP = PS I 
3 1 0 C X = 0. 

Y = P I HA, I F 

CALL SEARCH! LAMBDA, PS I , X , Y , D EGF ) 

IF (KILL.E0.2) GC TC 3200 

CALL SHIFT! SAVFL, SAVEP .LAMBDA, PS I , 1 , DEGF , D I ST ) 

GO TO 1 6C0 
3200 SA VEL = L AMBDA 
S AV EP= PS I 
GO TO 1600 

3300 DEGLAM=LAMBDA*RADDEG 

IF (DEGLAM.GT .180. ) GO TO 3400 
X = SYMBOL ! 3 ) 

GO TO 3500 

3400 DEGLAM=360 . -DEGLAM 
X= SYMBOL ( 4 ) 

3500 NDEGX=DEGL AM 
FDEG=NDEGX 

NMI NX=(DEGLAM-FCEG)*60.+.5 
DEGPHI=ATAN( FNCECC*TAN( PSI ) )*RADDEG 
IF (DEGPHI-LT.O. ) GC TO 3600 
Y= SYMBOL ( 1 ) 

GO TO 3700 
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3600 Y= SYMBOL (2 ) 

DFGPHI =— DEGPHI 
3700 NDEGY=DFGPHT 
FDFG =NDFGY 

N M I N Y = ( DEGPHI-FDEG )*60.+. 5 

WRITE (6,2) NO FGX , N MI NX , X, MDEGY , NM IN Y , Y 

2 FORMAT ( • O' ,30X, 13 , '-' , 12, A1 ,' , ' , 13 , •-' , 12 , Al, • IS THE 
1' LF AST SQUARES SOLUTION.') 

WRITE (6,3) OEGF 

3 FORMAT ( 31 X , ' OBJECTIVE FUNCTION =«,F8.2,' DEGREES', 

1' SQUARED.') 

WRITF (6,4) NUMBER 

4 FORMAT ( 3 0 X , I 3 , ' DEGREES OF FREEDOM') 

WRITE (6,5) 

5 FORMAT {'O' ,45 X, ' SOLUTION SIGNAL SITE',/,30X, 

I* ID INTERCEPT AZIMUTH AZIMUTH RANGE PLOT', 

2' THETA VALS ' , / , 30 X , ' NR SITE NAME BETA - THFTA ', 
3'= OMEGA NM NORM NM AZ • ) 

S U M A R C = 0 . 

DC 4000 K = 1, NUMBER 
I D= SI TF ( 7 , K ) + . 5 
DTHET A = THETA ( K ) *RADDEG 

CALL AZD I ST (LAMBDA, PS I , S IT E ( 1,K ) ,SITE(4,K) , SITEAZ, 

IDT ST) 

NMSITE=DIST+.5 
ARCC=(DI ST /3 4 39. 92) 

ANGL EB=ABS ( DCMEGA ( K ) )*DEGRAD 
ARCB=ARS IN (SIN ( ANGLES )*S INC ARCC ) ) 

ANGLE A=A RCOS ( T AN ( ARCS) *COTAN( ARCC) ) 

STAT=ARCC*RADDEG 

IF (ST AT. GT. 90 . ) ST AT= 180 ,-STAT 

IF (DOMEGA(K).GT.O. ) GO TO 3800 

AZTHET= SITEAZ- ANGLE A*RADDEG 

IF (AZTHET.LT.O.) A ZT HET=AZT HET+360 . 

STAAZ ( K ) =AZTHE T+90. 

IF (STAAZ(K) .GT .360. ) STAAZ ( K) =STAAZ ( K) -360. 

GO TO 3900 

3800 AZTHET=SITEAZ+ANGLEA*RADDEG 

IF ( AZTHET .GT .360. ) AZTHET=AZTHET— 360. 

STAAZ(K) = A ZTHE T-90 . 

IF (STAAZ(K).LT.O. ) STAAZ ( K) =STAAZ ( KH-360 . 

3900 NMTHET=ARCB*3439.92 

WRITF (6,6) ID,SITE(8,K) ,SITE(9,K) , S IT E ( 10 , K ) , BET A ( K ) , 
1 DT HET A , DOMFGA ( K ) , NM S I TE , AZ THE T , NMTHE T , ST A A Z ( K 1 

6 FORMAT (30X, 13 , 1X,3A4, F8.3, F8.3, F7.3, 16, F6. 1, 14, F6. 1) 
4000 SUM ARC = SUMARC+STAT 

DNUM=NUMBER 
AVGARC= SU MARC/ DNUM 
WRITE (6,7) AVGARC 

7 FORMAT ( 'O' ,30X, 'MEAN ARC STATISTIC =',F5.1,' DEGREES. 
COUNT = 0 . 

SUMANG = 0. 

N END=NUMRER-1 
DO 4100 1=1, NEND 
JB EG IN =1 + 1 

DO 4100 J = JBEGIN, NUMBER 
ACUTE =ABS( STAAZ ( I ) -STAAZ (J ) ) 

IF ( ACUTE .GT. 1 80. ) ACUTE = ACUTE- 1 80 . 

IF ( ACUTE. GT .90. ) ACUT E= 180 .-ACUTE 
SUMANG = SUMANG+ ACUT E 
4100 COUNT= COUNT+ 1 . 

AVGANG = SUM ANG/ COUNT 
WRITE (6,8) AVGANG 

8 FORMAT ( 3 1 X , 'MEAN ACUTE ANGLE OF INTERSECTION =',F5.1, 
1* DFGRFFS.' ) 

CALL PLOT (LAMBDA, PSI,DEGF, NUMBER) 

GO TO 0900 
4200 CONTINUE 
STOP 
END 
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SUBROUTINE SF ARCH ( L , P , X , Y , SA VE F ) 

C GOLDEN SECTION SEARCH 
REAL LAMBDA! 2), L 

DIMENSION DFL( 2) , FDEL ( 2) ,PSI ( 2) 

COMMON/SITES/S I T F ( 10,19) , NUMBER 
COMMON/ SIGNAL/ THETA (19 ) 

DATA T WOPI/6. 2 83 185/»PI HALF/1. 570796/, PI/3. 14 1593/, 

1 RADDEG/57. 29578/ 

SA VFL=L 
S AV EP= P 
I F L AG=0 
BEGIN =— 1 . 

END= 1 . 

J = 1 

1000 DEL ( J ) = 0 . 6 1 8*6 EG I N+ 0. 382* END 
GO TO 1200 
1100 I F L AG =2 

DEL ( J ) = 0 .382*8 EG I N+ 0. 6 1 8* END 
1200 E DEL ( J ) =0 . 

LAMBDA ( J ) -L + X* DE L ( J ) 

PSI< J) = P+Y*DEL ( J ) 

C FIND CORRESPONDING TEST LONGITUDE IN THE INTERVAL 
C (0.TWOPT). 

IF (LAMBDA(J) .GT.TWOPI ) LAMBDA ( J ) =LAMBDA { J J-TWCP I 
IF (LAMBCA(J).LT.O. ) LAMBDA (J)=TWOPI+LAMBDA(J ) 

C FIND CORRESPONDING TEST LATITUDE IN THE INTERVAL 
C ( -PI HALF, PI HALF ) . 

IF (PSI ( J) .GT. PIHALF) PS I ( J ) =P I-PS 1 ( J ) 

IF (PSI ( J ) .LT. -PI HALF) PSI ( J) =— PI — PSI (J) 

COS LA M = CCS ( L AM BC A ( J ) ) 

SI NLAM=SI N ( LAMBDA! J ) ) 

CCSPSI=COS (PSI ( J ) ) 

SI NPS I =S IN ( PS I l J ) ) 

VARLAM=LAMBDA( J) 

DO 1700 K = 1, NUMBER 

CALL AZCCMP (COSLAM, SI NLAM, COSPS I , S INPS I , K , VARLAM, 

1 A , B , C ) 

TFST=A*'B4C-THETA(K ) 

IF (ABS(TEST)-PI ) 1600,1600,1300 
1300 IF (TEST) 1400,1600,1500 
1400 C = C+T WOP I 
GO TO 1600 
150 0 C=C— T WOP I 
1600 CMEGA=A*B+C-THETA(K) 

1700 FDEL( J) =FDEL(J)+CMEGA**2 

IF ( IFLAG.GT.l ) GO TO 1800 
J=2 

GO TO 1100 
1800 CONTINUE 

C RETAIN VALUES CORRESPONDING TO CURRENT MINIMUM, ADJUST 
C INTERVAL AP'D TEST NEXT DELTA SYMMETRIC WITH CURRENT DELTA. 

IF ( FDEL ( 1 )-FDEL ( 2 ) ) 1900,1900,2000 
1900 END=DEL ( 2) 

PS I ( 2 ) = PSI (1) 

LAMBDA (2 ) = LAMBCA(1 ) 

FDEL(2)=FDEL(1 ) 

DEL ( 2 ) =D EL ( 1 ) 

J = 1 

GO TO 2100 
2000 BEGIN=DEL ( 1 ) 

PSI ( I ) =PSI (2) 

LAMBDA ( l)=LAMBDA(2) 

F DEL ( 1 1 =FDEL ( 2 ) 

DE L ( 1 ) =DEL ( 2 ) 

J = 2 

C TEST SIZE OF THE INTERVAL 
2100 IF ( IEND-BEG1N) .LE.0.01 ) GO TO 2200 
GO TO ( 1000, 11 00 ) , J 
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2 200 



MIN=1 

IF (J.F0.1) Ml N = 2 
L=LAMBPA(MIN) 

P=PSI ( MIN) 

DELTA=DEL( MI N) 

DFGF=FOELC M I N ) * P. ADDEG** 2 
IF (DEGF.L. E.SAVEF) GO TO 2300 
L = SAVE L 
P = SA VEP 
RETURN 
2300 SAVEF=DEGF 
RF TURN 
FND 



SUBROUT I NE A ZD I ST ( RL AM , RPS 1 , TLAM »T PS I , A Z , D I ST ) 

C AZDIST COMPUTES AZIMUTH AND DISTANCE OF A TARGET WITH 
C RESPECT TC A REFERENCE POINT 

DOUBLE PRECISION C 0 SR P , SI NRP , COST P , S I NT P , UC , RNUM , UG , 

1 DFN < COS B , R EFLAM » RE FPS I « TAP LAM , TA RP SI , AR C 
DATA TWOPI /6. 2 83 185/, R ADDEG/ 57 <29578/, P 1/3 , 141593/ 
REFLAM=RLAM 
REFPS 1 =RPS I 
TARLAM=TLAM 
TARPSI =TPSI 
COSRP=DCOS (REFPSI ) 

S I NRP =DS I N ( REF PS I ) 

COSTP=DCOS ( TARPSI ) 

S I NT P = D S IN (TARPSI ) 

UC-DCOS { TARLAMJvCCCSI REFLAM )+DS INI TARLAM )*CS I N( RE FLAM) 

RNUM=S 1 NTP-COSRP--CO STP- : S I NRP^UC 

UG=SINTP*SINRP+COSTP*COSRP*UC 

IF (UG.GE. 1.0D3) UG=1<0D0 

IF (UG.LE.-l <ODO) UG=-1.0D0 

ARC= CARCCS t UG ) 

DEN=DS 1 N ( ARC ) 

IF (DABS (ARC ) -LT.O. ID— 4) GO TO 1600 
IF (CABSCRNUM) .GT.O.OD-IO) GO TO 1000 
COSB = 0- ODO 
GO TO 1100 
1000 COSB=RNUM/DEN 
1100 CONTINUE 

CCSB=RNUM/ DEN 

IF (COSR.GE. 1. ODO) C0SB=1.9D0 
IF (COSB.LE.-l. ODO) COSB=-1.0DO 
BFTA=CARCOS(COSB) 

TF ST =REF LAM-TA RLAM 
IF (TEST) 1200,1500,1300 
1200 IF (TEST+P1) 1400,1500,1500 
1300 IF (TEST-PI) 1400,1500,1500 
1400 BET A=TWOP I -BET A 
1500 A Z = BET A^RACDEG 

DI ST=ARC"3439. 92D0 
GO TO 1700 
1600 AZ = 0 . 

DI ST=0. 

1700 CONTINUE 
RETURN 
END 



SUBROUTINE A ZCOM P ( CCSL AM , S I NL AM , COS P S I , S I NPS I , K , 

C AZCOMP COMPUTES POINT AZIMUTHS WITH RESPECT TC SITES 
1 LAMBDA , A , B , C ) 

REAL N, LAMBDA 

COMMON/SITES/SI TE( 10,19) , NUMBER 

DATA P 1/3. 141593/ ,T WOP 1/6.28318 5/, SCREEN/ .5599999/ 
A = 1 . 

C=0. 

UC=COSLAM*SITE (3,K ) +S I NL AM*S IT E ( 2 , K ) 

N = S1 NP S I * S I T F ( 6 , K ) -COS P S I* S I T E ( 5 , K ) *UC 
UG=S INPSI*SITE ( 5,K )+C0 SP S I* S I TE ( 6 , K ) *UC 
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IF (UG.GE.l.) UG=SC RE EN 
IF (UG.LF.-l.) UG=- SCREEN 
D=S IN ( ARGOS ( UG ) ) 

CGSB=N/D 

IF (COSR.GE.l.) CO S B = SC RE E N 
IF (C0S8.LF.-1. ) CO SB=- SCR E E N 
B-ARCC'S ( COSB) 

TE ST = S I TE ( 1 ,K)-LAMBDA 
IF (TEST) 1000,1300,1100 
1000 IF (TEST+PI) 1200,1300,1300 
1100 IF (TFST-PI) 1200,1300,1300 
1200 A=— 1 

C=T WOP I 
1300 RETURN 
END 



SUBROUTINE SHI FT (FPOML , F ROMP , TOL , TCP , NRTYPE , DEGF , DI ST) 
C THIS SUBROUTINE COMPUTES AND WRITES SEARCH PROGRESS 
DIMENSION T YPE S ( 2 ) 

DATA RADDEG/57. 29578/, FNC ECC / 1 , 006768/ 

DATA TYPES/' GOLD’ ,' GRAD'/ 

CALL AZDI ST(FROML,FRCMP,TOL , TOP , AZGO , DI ST ) 

D E GL AM= TOL *R ADD E G 

DFGPHl =ATAN( FNCECCYTAN (TOP) )*RADDEG 
IF ( NRTYPE .EO. 2) GO TO 1000 
TYPE=TY P ES ( 1 ) 

GO TO 1100 
1000 TYPE=TYPES( 2 ) 

1100 WRITE (6,1) TYPE, D I ST , AZGO, DEGLAM, DEGPHI , DEGF 
1 FORMAT (31X, A4 ,F7.1 , ' NM ON AZ',F6.1,' TO » , F8 . 3 , ' E • , 

1F8.3, ' . F=« ,F10. 3) 

FROML=TOL 
FROM P = T OP 
RETURN 
END 



SUBROUTINE Pl.OT( LAMBDA, PSI , DEGF , NUMB ER ) 

C PLCT COMPUTES AND WRITES THE CONFIDENCE REGION 
REAL LAMBDA, L 

DOUBLE PRECISION VAR PS I , P I HALF , P I , ARC , COS COL , COL AT , 

1 CGSANG 

COMMCN/SIGNAL/ THETA (19 ) 

DIMENSION OUT ( 52) , SYMBOL ( 3) 

DIMENSION CH 19 00 (19 ) 

DIMENSION JSCALE(ll) 

DATA SYMBOL/'-' , ' • , 'S'/ 

CAT A PI HALF/ 1 . 57 0796326 DO/ , P I / 3 . 1 4 1 59265 3D0/ 

DATA DEGRAD/1 . 74 53 2 9E-2/ , T WOP I /6. 2831 85/ , 
1RA0DEG/57.29578/ 

DATA CHI 900/2. 71, 4. 61, 6. 25,7. 78, 9.24, 10.6, 12.0, 13.4, 
114.7,16. ,17.3, 13.5,19. 8,21.1,22.3,23.5,24.8,26. ,27.2/ 
BCUND=CHI900 (NUMBER ) 

PROB =. 900 
NS I ZF = 400 

0900 S I 7.E=NS I ZE/2 

T0P=PSI+1. 05vSIZE*2.9070443E-4 

KCUNT = 0 

MARKV=0 

DO 2700 I = 1,41 
ADD= 0 . 

IF (MARV.NE.2) GO TO 1100 
DO 1000 J = 1,51 
OUT ( J ) = SYMBOL ( 1) 

1000 CONTINUE 

GO TO 2500 
1100 T AK F= I 

DIRECTS . 

VARPSI=TOP-TAKE*0. 0 5* S I ZE*2. 9070443E-4 
IF (VARPSI .LT.PIHALF) GO TO 1200 
VARPSI =PI-VARPSI 
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D IRECT =— 1 . 

ADD = P I 
GO TO 1300 

1200 IF (VARPSI .GT.— PI HALF) GO TO 1300 
V ARPS I =— P I —VAR PS I 
D I RFC T=- 1 . 

A00= P I 

130 0 PARC=1 . 04 s ? S I ZE#2 .9070443E-4 
MARKH= 0 

00 2500 J = 1,51 

IF ( MARK H. EC .2 ) GO TO 2100 

T A K E = J 

ARC=PARC-TAKE*0.04-*SIZE*2.90 70443E-4 
COSCOL=DCCS( ARC) *DSIN( VARPS I ) 

IF (COSCOL.GT. 1 .DO) C.OSCOL = 1.DO 
IF (COSCCL .LT.-l .00 ) COSCOL=-1.DO 
COLA T = D A PC GS (COS COL) 

P= PI HALF-COL AT 

CCSAKIG = DCOT AN (VARPS I )*DCOTAN ( COL AT ) 

IF (COSANG. GT. 1.00) C0SANG=1.D0 
IF (COSANG. LT.-l .00 ) C0SANG=-1.D0 
ANGLE=DARCCS (CCSANG ) 

IF (J.LE.26) GO TO 1400 
L = L AM80A + D I RECT y ' ANGLE + ADD 
GO TO 1500 

1400 L = LAMRDA— 0 I RFC T“ANG LE+ ADD 
1500 CCNTINUF 

IF (L.GT.TWOPI ) L = L -T W 0 P I 
IF (L.LT.O.) L =T WO P I +L 
COS L AM = COS ( L ) 

SINLAM=SIN(L) 

COSPSI=COS ( P ) 

SINPSI=SIN(P) 

DEGF=Q. 

CO 2000 K = 1, NUMBER 

CALL A ZC CMP (COSLAM , SI NL AM , COSPS I , S IN PS I , K , L , A , B, C ) 
TFST = A*B-i-C-THETA<K) 

IF (ABS (TEST )— P I ) 1900,1900,1600 

1600 IF (TFST) 1700,1900,1800 
170 0 C = C + T WO P I 
GO TO 1900 
1800 C^C-TWOPi 
1900 OMFGA=A*B+C-THETA( K ) 

2000 CFGF=DEGF+(CMEGA4RArDEG)**2 
2100 OUT(J) =S YM80L ( 1 ) 

MARKE R=MARKH+1 

GOTO (2300,2200,2500), MARKER 
2200 IF { DEGF.GT. BOUND) MAR KH =2 
2300 IF ( DEG C .GT .80UN C ) GO TO 2400 
OUT ( J ) =SYM Q GL ( 2 ) 

K0UNT=K0UNT+1 
MARKH=1 
N A R K V = 1 

IF (I .GT. 1 ) GO TO 2350 
NSIZE=NS IZE + 20 0 
GO TO 0900 

2350 IF (J.GT.l) GO TO 2400 
NSIZE=NSIZE+NSIZE 
WRITE (6,7) 

7 FORMAT ( ’ ) • ) 

GO TO 0900 

2400 IF (I.NE.21) GO TO 2500 

IF (J.E0.26) OUT (J )=SYMB0L(3) 

2500 CONTINUE 

IF ( MARKV.NE. 1 ) GO TO 2600 
IF (MARKF.EO.O) M A R KV = 2 
2600 CONTINUE 

IF ( I .NE. 1 ) GO TO 2650 
WRITE (6,1) PROB,NSIZE,NSIZE 

1 FORMAT ( ‘O' ,30X, F6.3, ' CONFIDENCE REGION PLOT ON’, 15, 

1 • NM BY’ , I 5, ’ NM AREA' ) 

WRITE (6,3) 
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3 FORMAT ( 37X, 17 ( ’ » ) , * — » ) 

2650 CONTINUE 

VF RT=? 1 — I 

NM=A3S( 0.05* VERT*S 1ZE )+.5 

IF (I.GT.21) NM=-NM 

WRITE (6.4) NM, ( CUT ( J) ,J = 1,51) 

4 FORMAT (‘ ’ ♦ 30 X . I 5 . ' I'.51A1,'I') 

2700 CONTINUE 

WRITE (6.3) 

DO 2800 J = 1,11 
ADO=J 

JSCALE ( J) =ABS( (0.2*ADD-1.2K$IZE) + .5 
2800 IF (J.IT.6) JSC Al E ( J)=— JSCALE( J) 

WRITE (6,5) (JSCALE(J) , J = 1 , 1 1 ) 

5 FORMAT (34X.11I5) 

IJAREA=G.002 A S IZF**2+.5 
N MAR EA = K CUNT* I J A R E A 
WRITE (6,6) NM AREA 

6 FORMAT ( 30 X , ! CONFIDENCE REGION AREA =»,17, 

1» NAUTICAL MILES SQUARED.') 

RETURN 

END 
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GOLD 


4398.4 


NM 


ON 


AZ 


270.0 


TO 


286. 740E 


0.0 . F = 


20254.277 


GOLD 


4675 . 4 


NM 


ON 


AZ 


180.0 


TO 


2 86 . 740 E 


-77.953. F= 


11288.250 


GOLD 


1021.9 


NM 


ON 


AZ 


134. 1 


TO 


16. 323E 


-77.953. F= 


8562 .371 


GOLD 


1098 .3 


NM 


ON 


AZ 


0.0 


TO 


16.323E 


-46.449. F= 


3712.686 


GOLD 


1108.5 


NM 


ON 


AZ 


99.8 


TO 


43 . 155 E 


-46.449. F= 


868.700 


GOLD 


412.8 


NM 


ON 


AZ 


0. 0 


TO 


43. 155E 


-39.571. F= 


1 1 1 .312 


GRAD 


271 .3 


NM 


ON 


AZ 


94.0 


TO 


49.004F 


-39.741. F= 


15.526 


GRAD 


84. 6 


NM 


ON 


AZ 


53.8 


TO 


50.461 E 


-38 .898 . F= 


0 .158 


GRAD 


10. 2 


NM 


ON 


AZ 


76. 8 


TO 


50. 673E 


-33.859. F = 


0.050 


GRAD 


0.1 


NM 


ON 


AZ 


52.7 


TO 


50.675E 


-38.858. F= 


0.050 



50-40E, 38- 5 1 S IS THE LEAST SQUARES SOLUTION. 
OBJECTIVE FUNCTION = 0.05 DEGREES SQUARED. 
3 DEGREES OF FREEDOM 



ID INTERCEPT 
NR SI TF NAME 
9 HONOLULU 
19 WICHITA 
10 LOS ANGELES 



SOLUTION SIGNAL 
AZIMUTH AZIMUTH 
BETA - THETA = 
228.098 228.241 
102.627 1 02 . 797 
121.457 121.424 



SITE 
RANGE 
OMEGA NM 
-0.142 9001 

-0.170 9301 

0.033 10199 



PLOT THETA VALS 
NORM NM AZ 
27.1 4 117.1 

187.3 4 277.3 

335.1 0 245.1 



MEAN ARC STATISTIC = 21.8 DEGREES. 

MEAN ACUTE ANGLE OF INTERSECTION = 34.7 DEGREES. 



0.900 CONFIDENCE REGION PLOT ON 400 NM BY 400 NM AREA 
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■- 5^ rt .v * £ 



^ t> 3V. a f' ^ ft .V ii 4“ 't: ft ft ft ft ft: 1 : ft ft *£■ r> r J 

^ r ft ti; r? -V. if ^ ^ >* ft ft % •’< * At !• / i ? . P. t: £ ft J 

ftr ft ft: ft ^ xk rt & ftzftft ftc’Zzft ft. ft-: ft: ft ft: :ft ft ft? ft' ft zk ft J 
ft: k ft* ft\ fti -f: *, *> A zk ft ft~: >t ftft A ^ $.&:):&. ft P J 
ftftft: ly*': & i\* :*: rCst ft *. k :?. ftt* #ftz A $ 4 i ft-. r v J 

:ki\zk^rk *f & t kfiftk? ft ft-*'. ftzftkzkft- frftk±*r.t > ;> i< J 

:*; t .v A ^ k y r> * >• ^ ^ ^ ^ ^ ^ t~ ft -ft ftr £ :• rVN^ 

p\ ft; &ftr ft: p it 7 >. ftx tk :ft ft 'M if . k jir ^ &&$■$. & $i ft; 
ftz * ftr. i. i h *: v> ft.ft&tfrt* tS A 4. 4 ft- :i tft a 4? tV ft ft-zk ft ft ftz ft J 

ijr^ st A £- k 1 : k rV il w jfe tV .t- fir * ^ :£r ir t a r £ i tV + J 

'Aft' A rk t> rt: ftr. :kft: r: ft? ft: tt r> tkft£ rirA. J 

45 f- j?. -A'i&ftr.x ft- ft; A A t ^ ft P A A ^ ^ 5 *. 4: ^ r: j 

ft, ft ft* :;• A'i;P ft: *:*; ft # ftr. f : . ft ft. ft. ft & ft ft is. ft. -i ft- ^ ^ i ft if J 

ft, ft: ft ft ft. A A A A ft ft ft. ft ?> ft ft ft. A ft ft <: ft ft A ft A ftftc J 
ft A A ft. ft ft ft: ft ftA t ■*?. ^ A ft & £ * A st ft. ft ft ft ft J 

ft rk A ft. ft ft ft ft t': z'uti ft ft ft A ft A ft ft ft ft t ft. 1 ft ft ft ft. ft J 
5^ ft: ft. ft: ft t ft. ft A. ft. ftfuAft Aft ft : ft A ft ft ft ft 'ft. ft ft k ft Aft J 



-200 -160 -120 -80 -40 

CONFIDENCE REGION AREA = 



0 40 SO 120 160 200 

6080 NAUTICAL MILES SQUARED. 



97 



GOLD 


6104.5 


NM 


ON 


AZ 


270. 0 


TO 


258. 323E 


0.0 


F = 


86.016 


GOl D 


276.9 


NM 


ON 


AZ 


180.0 


TO 


258.323E 


-4.643. 


F = 


79.381 


GRAD 


2293.2 


NM 


ON 


AZ 


167. 1 


TO 


268.963 E 


-41 .796 . 


F= 


10 . 165 


GRAD 


1199.7 


NM 


ON 


AZ 


160. 2 


TO 


282. 3 OOE 


-59.978. 


F = 


0.043 


GRAD 


77.3 


NM 


ON 


AZ 


151.7 


TO 


283. 557E 


-61.102. 


F = 


0.013 


GRAD 


0. 9 


NM 


ON 


AZ 


97.3 


TO 


283. 589E 


-61.104 . 


F = 


0.013 


GRAD 


0.0 


NM 


ON 


AZ 


0.0 


TO 


283. 590E 


-61.104. 


F = 


0.013 



7 6-2 5W» 61- 6S IS THE LEAST SQUARES SOLUTION. 
OBJECTIVE FUNCTION = O.Ol DEGREES SQUARED. 
3 DEGREES OF FREEDOM 



ID 


INTERCEPT 


SOLUTI ON 
AZIMUTH 


SIGNAL 
AZI MUTH 




S ITE 
RANGE 


PLOT 


THETA 


VALS 


NR 


SITE NAME 


BETA - 


THETA = 


OMEGA 


NM 


NORM 


NM 


AZ 


7 


EL PASO 


165. 841 


165.922 


-0.081 


5747 


244 .6 


4 


334.6 


17 


RENO 


159.655 


159. 655 


0.000 


63 8 3 


56. 4 


0 


326.4 


8 


HELENA 


162.349 


162.270 


0.080 


6673 


64. 5 


4 


334.5 


MEAN ARC STATI STIC = 75 
MEAN ACUTE ANGLE OF INT 


.6 DEGREES. 
ERSECTION = 5. 


5 DEGREES. 







0.900 CONFIDENCE REGION PLOT ON 2400 NM BY 2400 NM AREA 



1200 
1140 
1080 
1020 
960 
900 
840 
780 
720 
6 60 
600 
540 
480 
420 
360 
300 
240 
180 
120 
60 
0 

-60 
-120 
- 180 
-240 
-300 
-360 
-420 
-480 
- 5 AO 
-600 
-660 
-720 
-780 
-840 
-900 
-960 
-1020 
-1080 
-1140 
-1200 



I it t ;'r • iyjt \ ^ s'; - ' r :•£ “u S' ; f.. A- £• * >; i: ft. J? ij: >•: i~. A J 
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Jr? 
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I :-. t ;*: yt rf jL; n; A* S. vc r» rf\ ^ rf* j v A rf rf A rf A rf rf. A A rf rf 

I rf.: J :rf*:5;:.r.:f:rf-rf ^r v .^:^rf2?fc^74:rf.*AA A A rfi Arfrrf^rf 
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T rf. ^ rf rf; ‘u rfrf rf n*:|: v A rf: ^ ; * rf rf rf- rf rf rf rf rf A **\7' a 

J ^ *• rf rf f tt y A rf rf sV rf rf rf. rf rf rf*-'* rfrf rfrf, r> rf r<: ‘.A: rf: A : i: u rf r?: 7*^ 
l & % ff. :Vrf rkrfrf rf rf rf A i& t % rf *A sf, * rf ^ rf ft rf ^ f-; rf ^ :fc rf : 

X rf rf ■ :': rf vx izia A rf, rf rf; 15* rf rf r> rf:f. .** S:y.rfrf^rf rf rfrft;^rf>Ar^A^ 

X rf i; s; as rf .*£ £ rf f rf; rf rf ft £ rf rf* rf 5 V ft rf ft j^r ft rf rf: ft rf: rf ft rf rf ft rf r$t ^ 

X : « * :s ft ft ^ ft ^~A' r: ft ft :< ft ft ft }’.: ft ^ :l- ft ft ft rf 2/ ft 'h ft ft r** ft ft ft ft f ft ft 

X 5 >^ftft 5.^ rfrfrfft ^ Aft rfrf^rfrfftrfftc^ i 



Ly r rf cSi 2f;rf rl*. rf.ft V rf ft ft ft; ft; A ?V’ rf •*» A k k ft X 
rf :‘5 ft ft: rf • ;; rf K? ft : jJ % rf rf A’ rf ft; rf ft: ^ $ 5f rf I 
rfz 'ft -rfrfrfft rf *: tc rf *♦ rf rfst ft;rf.ft. X 
rf rf ft rf ftr ft; A A ft ■$, rf * rf A rf :< ft. r ' rf y. ft rf ft' T 
rf rf rf if t \k y? y rf. ft: rf •': ^ rf ♦ rf rf ft rf :>ft- J 
:' rf **.rf rf >v' j-hrf ft - rf 7^' rf rf ft. rf.rf rf rf rv*: V X 



ft ;•< rf c rf t rf «?- ft rf ir rf rf ft ft k: ri i X 
rfftc ft; ft: ft * 4. ft rf ft: ft *f ft: ft*, h ft- rf* rf y : X 
rt yz rf rf :r*.:^ ft ft- ft. rf :V v A jp aft rf rf A J 
rf. A rf :rf rf £ rf rf ^ A rf V A X 
Axbr^^-^rfrfrfft rf A X 
u < rf rf ^ £ rf rf ft: rf £ rf 5L. ft'r J 
y ft-, J: rf .v rf rf rf ft ft" A* rf rf rf ft* X 
rf :-s ft- : > i: rt z> rf 2? ft rf J 

ft vfc rf ^ A k rf rf : rf : ^ rf ^ f 

. 7 z S' A rf rf ft; ft: ft* rf rf X 
rf ik r -: ft* A ft: ft ft: ft T 
rf 4- rf - ' ft' A y .< rfft* f 

. ft rfrf rfrfrf X 
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-1200 -960 -720 -480 -240 0 240 480 720 960 1200 

CONFIDENCE REGION AREA = 311040 NAUTICAL MILES SQUARED. 



98 



GOLD 


7810 . 8 


NM 


ON 


AZ 


270.0 


TO 


229 .901E 


0 .0 . 


F = 


663.993 


GOLD 


995.3 


NM 


ON 


AZ 


180. 0 


TO 


229. 901 E 


-16.684. 


F = 


13.884 


GOLD 


97 .0 


NM 


ON 


AZ 


90.2 


TO 


231. 588E 


-16.684. 


F = 


3.408 


GOLD 


3 . 0 


NM 


ON 


AZ 


0.0 


TO 


231.588E 


-16.620 . 


F= 


3.381 


GRAD 


10.1 


NM 


ON 


AZ 


67.7 


TO 


231. 750E 


-16.556. 


F = 


3 .289 


GRAD 


0 .0 


NM 


ON 


AZ 


0.0 


TO 


23 1 . 7 50 E 


-16.556. 


F = 


3.289 



128-15W, 1 6-3 3 S IS THE LEAST SQUARES SOLUTION. 
OBJECTIVE FUNCTION = 3.29 DEGREES SQUARED. 

4 DEGREES OF FREEDOM 



ID 


INTERCEPT 


SOLUTION 

AZIMUTH 


SIGNAL 

AZIMUTH 




SITE 

RANGE 


PLOT 


THETA 


VALS 


NR 


SITE NAME 


BETA - 


THETA = 


OMEGA 


NM 


NORM 


NM 


AZ 


14 


NASHVI LLE 


224.264 


224. 377 


-0. 113 


3929 


306.1 


6 


36.1 


2 


ATTU 


124.892 


126.340 


-1.448 


5135 


238.9 


86 


328.9 


9 


HONOLULU 


140.026 


138.983 


1 .043 


2854 


50.6 


46 


320.6 


16 


NORFOLK 


232.077 


231.774 


0.303 


4406 


130.0 


17 


40.0 


MEAN ARC STATISTIC = 68.0 DEGREES. 

MFAN ACUTE ANGLE OF INTERSECTION = 50.9 


DEGREES . 







0.900 CONFIDENCE REGION PLOT ON 400 NM BY 400 NM AREA 



200 I £ r; A k A A a £ :V- si- 
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10 I J > ^ <* ^ ^ # -•• A A \ r' jC: :** A v 

0 I V ^ t ^ ^ •?' 1 1 *• > ^ JV* * ifc A 

— 1 0 J A :•). A r'- :*? A :>• A ^V: ^ ^ 

20 j !«*• *5i t*' , A t* t f r A: A A i" '*• 

— 30 ] V i z A ^ iv :•: A i : :y * <; :' 

— 40 I % ^ ^ ■ A ^ ^ A 

— 50 I ^ ^ rjr 

— 6 0 t ^ rt :?; :l::*r A *i: A f: iV st tV : 

— 70 J **: A r> 'Jzr%:tk': :\'W <c 5z ?- v 

— 8 0 I<*£ t* * J ^ * r? T *>A: A •*?. :> A ?'< ^ r/ 

— 9 0 I :J; r?t * # *' r y‘f A ••" v* .> t A A v 

-100 1 :”r A rV y‘ >: r/% > rA > * ^5 j?: A r^ ? * 4 . s/ ,*. ; . 

"110 1 ^ V -V- w ?[.: A A r»:!> A ^ ^ 

" 120 I ^ V V-! ^ A*Vrt^; j 

"130 I rv;?*:* ,v H:r72 :^::*< 

-140 I A ** u v :’t!.> ^ A t * i 

“150 I ’ •/ ^ J ^ :> 2 f* A rye v\k A A ;> 

"160 1 ^ fc 

-170 I v L jf ' 12 ’ r ^ 1 A I T >*4 ^ / '/^ 1 ^ y jl? 3f 

"180 I **-' ^ •t’ * * A A. 

" 190 I ^ y^r-::k^ Ar? 

"200 I <2 ^ S' rk tiyy'rC 4 * r\- A i > v ! . 4 



r A . r* A ^ A 4; A iA A A A rf: ?S? A ?< *' A A* iyC: 5c iji A ^ t /. r’ J 

r- ’> # ^ A ^ A %t <s ?> : -i> rfe:*- * 3^. A ;V t;< A A A- ij ^ < J 

# *T? * *ir. ;*. .^f: A !k ?V A r v . A *A .<? A A .'<’ w 4 . A d? A rAr i? r A A J 

jv *! x & jfr- y> •*.- ife ««: <2 i: :ir d. s> «£ ^ 47 ^ $ r> ^ ir ^ •*’ *• /c : rV A s'-. A J 

A ‘‘ rtS; : Art* AA A i» As’.iiniy ^ A J 

v. ; f ; ?> A A A A £ * A 1 i'z A jC ^ A**' l ; ’ A A?i £. :y^ i). ^ A fr fr ; J 

** AA AAfr $ fr-fr.A:Afr.y;:^: A . ^.. A* >, ^ ] 

:;V:A?VA :VAA Afr iA ^* I 

fr IfrV^ fr fr Jiff fr* *?* ^ :V 3jS fr fr <: 3ri i\ fr/N fr' >V fr. fr/^ <5 fr fr ^ ^ fr Zk X \ 

• fr Afr fr k fr- A Afr it m fr :A? fr frn fr A<> «*t fr ’ :V % te iy fr fr A. J 

it ^ fr fr A *>. fr fr 'A 4 ;f fc : v :y4 A :fj r : J :J/£' A fr r? A it* fr fr fr fr . fr fr , J 

£ r> fr fr & ?yi fr & r fr/-; fr fr ^ # > fr.. i; a. 4 s? c- j 

*4 :Vv ; A fr ?k a/:-, fr fr : •*. | :>::^fr fr-.d fr'frfr j 

i fr: rf^fr frfr . fr fr rx A A •: & r|: fr ^ fr/ “ \ 

-A ±:A * Sz fr •/ fr. fr: 4. 1/ :l A ;* ,A- fr fr fr J 

: *: Ar fr- A fr -fr: A ?■ i. ^ fr J fr-fr : . fr-fr frx J 

fr * & $ ■?. fr ^ . fr 1- $ 4 ;V fr v Sv j ? * J 

fr fr 4=. vi. fr fr A fc *; j£ ■? ^ fr* fr frfr.fr] 
:k A fr fr' 'A fr Afr fr: fr :\ fr fr fr fr fr li: ] 
► *4. } r •; A t i‘2 < fr’> I 

fr :y- fr s; fr fr fr. £ fr r ; fr fr fr fr t w r w . ] 
fr frfr frfrfrfrfr frfrnJ fr^-'frfrs, > ] 
frfrfrfr A Afr^frfrfrA&frfr fr ] 
raV: 2* fr fr fr fr fr tit fr fr > fr fr fr fr *' ] 
*:fr:tfrfrfrfrA> J 
fr fr A fr\: fr fr •». £ fr fr A fr A fr fr A A I 
frA.AfrfrNrAfr Afr.fr AfrA:; fr :• :’ T 
frfr^Afrfr ^ ^ i 

fr * fr:-t fr fr fr~ fr frH^r ❖ fr • A ?5 fr fr ^T* ] 

Vfr frA^frs>:Afr frfr\.^AA^frfr frrrA ] 
5-j fr; fr A A fr fr* fr fr. :k:U± fr fr ^ fr fr :*c fr fr. frVrf A A * :J; ] 

frx fr: :\ fr-j A fr fr fr; ± :- A ir A fr* frr A rJi fr. frfr fr frei* ft A: H fr Afr* ] 

A fr \ A b'i fri :\ A 4: fr' ly A A fri/i fr:A fr- fr. fr fr r i ^ A\fr- K A fre j fr ] 

fr* fr: fr. 4i* A frs A fr7JV A frs frrfr frk ?•- A f; j.*r i( fr fr: fr: fre A ] 

fr fr < fr. ^ V: fr >■? %. A i*:fr fr*fr :1* A r^* A fr y\ A . ywfr/AfrtA ] 

fr fr n fr A .V: fr* fr fr fr fr fr x \\ * *z fr fr ■* /’ ^ fr <7 < *: <• A r\?4 ;f*fr ] 

•& i' : A fr fr fr A AAA fr fr fr fr fr fr fr fr fr rfe fr • fr frXfr fr»* ] 
y -4- sSc s% A. fr 'frV fr A fr ;?* fr ?> fr fr fr : : \ frfr A A fr fr s fr 4: v T 

3* fr fr fr fr sfV: fr fr fr fr.fr: fr « fr xfrfrfr .y fr fr & fr frfr# fr: ^\ ] 



: r- fr fr:’: A fr fr * V^. :«: fr :’ A ?' . >zfr fr: fr fr. fr fr $ : 47 ♦* fr* fr: fr fr: fr. ♦ 

fr frfrrfr; fr’ fr v * s^Xs > A fr fr Afr ife fr £ £ fr fr: fr £ fr fr fr fr fr fr fr fr ] 



-200 -160 -120 -80 -40 0 40 80 120 160 200 

CONFIDENCE REGION AREA = 20880 NAUTICAL MILES SQUARED. 



99 



GOLD 


4399 . 1 


NM 


ON 


AZ 


90.0 


TO 


73. 272E 


0.0 . F = 


767.351 


GOLD 


1765.6 


NM 


ON 


AZ 


0.0 


TO 


73.272E 


29.574. F= 


26.527 


GOLD 


123.9 


NM 


ON 


AZ 


89.4 


TO 


75. 6 4 2 E 


29.574. F= 


4.538 


GOLD 


71 .1 


NM 


ON 


AZ 


0.0 


TO 


75. 642E 


30.763. F= 


3.487 


GRAD 


14.2 


NM 


ON 


AZ 


55.7 


TO 


75. 869 E 


30.897 . F= 


3.303 


GRAD 


0.2 


NM 


ON 


AZ 


40.6 


TO 


75. 8 7 1 E 


30.899. F = 


3.303 


GRAD 


0.0 


NM 


ON 


AZ 


0.0 


TO 


75.871E 


30.899. F= 


3.303 



75-52F, 30-54N IS THE LEAST SQUARES SOLUTION. 
OBJECTIVE FUNCTION = 3.30 DEGREES SQUARED. 
5 DEGREES OF FREEDOM 







SOLUTION 


SIGNAL 




SITE 








I D 


INTERCEPT 


AZIMUTH 


AZIMUTH 




RANGE 


PLOT 


THETA 


VALS 


NR 


SITE NAME 


BETA - 


THETA = 


OMEGA 


NM 


NORM 


NM 


AZ 


17 


RENO 


345.832 


345.356 


0.47 6 


6511 


102.9 


27 


12.9 


5 


DENVER 


359.185 


358.897 


0.288 


6589 


90.8 


16 


0.8 


12 


MIDWAY 


304.587 


304.956 


-0.369 


5329 


327 .7 


22 


57.7 


19 


WI CHIT A 
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MEAN ARC STATISTIC = 72.3 DEGREES. 

MEAN ACUTE ANGLE OF INTERSECTION = 28.0 DEGREES. 



0.900 CONFIDENCE REGION PLOT ON 400 NM BY 400 NM AREA 
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34-29E, 16-29N IS THE LEAST SQUARES SOLUTION. 
OBJECTIVE FUNCTION = 0.26 DEGREES SQUARED. 
5 DEGREES OF FREEDOM 
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MEAN ARC STATISTIC = 72.1 DEGREES. 

MEAN ACUTE ANGLE OF INTERSECTION = 28.4 DEGREES. 



0.900 CONFIDENCE REGION PLOT ON 600 NM BY 600 NM AREA 
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20-44W. 9-5 2 S IS THE LEAST SQUARES SOLUTION. 
OBJECTIVE FUNCTION = 3. 05 DEGREES SQUARED. 
6 DEGREES OF FREEDOM 
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MFAN ARC STATISTIC = 58.6 DEGREES. 

MEAN ACUTE ANGLE OF INTERSECTION = 22.0 DEGREES. 

0.900 CONFIDENCE REGION PLOT ON 400 NM BY 400 NM AREA 
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CONFIDENCE REGION AREA = 26720 NAUTICAL MILES SQUARED. 



102 



GOLD 


7158.7 


NM 


ON 


AZ 


90.0 


TO 


119.236E 


0.0 . 


F = 


318.482 


GO L D 


251.0 


NM 


ON 


AZ 


0. 0 


TO 


119.236E 


4.209. 


F = 


261.669 


GRAD 


1861 . 3 


NM 


ON 


AZ 


29. 7 


TO 


136. 486E 


30. 744. 


F = 


34. 133 


GRAD 


525. 6 


NM 


ON 


AZ 


47.2 


TO 


144. 45 IE 


36.464. 


F= 


8.503 


GRAD 


104.9 


NM 


ON 


AZ 


34.4 


TO 


145. 701E 


37.901. 


F = 


8.042 


GRAD 


5.3 


NM 


ON 


AZ 


222.2 


TO 


145. 627E 


37. 836. 


F = 


8.041 


GP AD 


0.6 


NM 


ON 


AZ 


222.0 


TO 


145. 61 8 E 


37.829 . 


F= 


8.041 


GRAD 


0.4 


NM 


ON 


AZ 


221.9 


TO 


145. 612E 


37.824. 


F = 


8.041 


GRAD 


0. 1 


NM. 


ON 


AZ 


221.8 


TO 


145. 6 1 1 E 


37.822. 


F = 


8.041 


GRAD 


0.0 


NM 


ON 


AZ 


0. 0 


TO 


145.611E 


37.822 . 


F= 


8.041 



145-37E, 37-49N IS THE LEAST SQUARES SOLUTION. 

OBJECTIVE FUNCTION = 8.04 DEGRFES SQUARED. 

6 DEGREES OF FREEDOM 







SOLUTION 


SIGNAL 




SITE 








ID 


INTERCEPT 


AZI MUTH 


AZIMUTH 




RANGE 


PLOT 


THETA 


VALS 


NR 


SITE NAME 


BETA - 


THETA = 


OMEGA 


NM 


NORM 


NM 


AZ 


16 


NORFOLK 


327.990 


328.189 


-0.199 


5685 


301.6 


11 


31.6 


4 


CORPUS CHRI S 


. 315.081 


314.619 


0.462 


5524 


142.2 


27 


52.2 


11 


MI AMI 


324.266 


322.103 


2.164 


6215 


132.2 


126 


42.2 


18 


SEATTLF 


298.750 


300 .372 


-1 .622 


3873 


319.2 


87 


49.2 


6 


DETROI T 


323.507 


324. 197 


-0. 690 


5328 


303.9 


41 


33 .9 


13 


MINNEAPOLIS 


316.823 


316.840 


-0.016 


4923 


307.8 


0 


37.8 


MFAN ARC STATISTIC = 80 


.8 DEGREES. 











MF AN ACUTE ANGLE OF INTERSECTION = 10.2 DEGREES. 

0.900 CONFIDENCE REGION PLOT ON 600 NM BY 600 NM AREA 



300 

285 

270 

255 

240 

225 

210 

195 

180 

165 

150 

135 

120 

105 

90 

75 

60 

45 

30 

15 

0 

-15 

-30 

-45 

-60 

-75 

-90 

-105 

-120 

-135 

-150 

-165 

-180 

-195 

-210 

-225 

-240 

-255 

-270 

-285 

-300 
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0 60 120 
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180 240 300 

MILES SQUARED. 



103 



GOLD 


2341 .5 


NM 


ON 


AZ 


9 0.0 


TO 


39. OOOE 


0.0 . 


F- 


5868.273 


GOLD 


575. 9 


NM 


ON 


AZ 


18 0.0 


TO 


39. OOOE 


-9 .656 . 


F = 


5288.270 


GOLD 


486 . 3 


NM 


ON 


AZ 


90. 7 


TO 


47. 215E 


-9.656. 


F = 


4666.781 


GOLD 


531 . 3 


NM 


ON 


AZ 


180.0 


TO 


47. 2 15 E 


-18.558. 


F = 


3902.349 


GOLD 


406 . 6 


NM 


ON 


AZ 


91. 1 


TO 


54. 3 55 E 


-18.558. 


F = 


3082.319 


GOLD 


426 .4 


NM 


ON 


AZ 


180.0 


TO 


54.355E 


-25.695. 


F = 


2172.056 


GRAD 


577.5 


NM 


ON 


AZ 


138.9 


TO 


61. 849 E 


-32.774. 


F = 


215.962 


GRAD 


95 . 2 


NM 


ON 


AZ 


36. 1 


TO 


62. 942 E 


-31.483. 


F = 


4.3 74 


GRAD 


10 .9 


NM 


ON 


AZ 


99.8 


TO 


63. 151E 


-31.514. 


F= 


2.152 


GRAD 


0. 1 


NM 


ON 


AZ 


25.4 


TO 


63.152E 


-31 .512 . 


F = 


2.152 


GRAD 


0.0 


NM 


ON 


AZ 


0. 0 


TO 


63. 152E 


-31.512. 


F = 


2.152 



63- 9E. 3 1-3 1 S IS THE LEAST SQUARES SOLUTION. 
OBJECTIVE FUNCTION = 2.15 DEGREES SQUARED. 
7 DEGREES OF FREEDOM 



ID 


INTERCEPT 


SOLUTION 

AZIMUTH 


SIGNAL 

AZIMUTH 




SITE 

RANGE 


PLOT 


THETA 


i VALS 


NR 


SITE NAME 


BETA - 


THETA = 


OMEGA 


NM 


NORM 


NM 


AZ 


5 


DENVER 


53.228 


53 . 230 


-0.002 


10049 


223.7 


0 


313.7 


15 


NEW ORLEANS 


100.649 


101.286 


-0.637 


942 5 


182. 5 


14 


272.5 


6 


DETROI T 


79.209 


78.237 


0.972 


9070 


32.3 


28 


302.3 


16 


NORFOLK 


90.386 


89.567 


0.820 


8779 


24.3 


27 


294.3 


7 


EL PASO 


91.177 


91.369 


-0. 192 


10276 


184. 1 


1 


274. 1 


17 


RENO 


342.343 


342 .395 


-0.051 


10306 


285.9 


0 


15.9 


8 


HELENA 


15.592 


15.892 


-0.300 


9875 


257.2 


4 


347.2 


MEAN ARC STATISTIC = 18.7 DEGREES. 

MEAN ACUTE ANGLE OF INTERSECTION = 42.9 


DEGREES . 







0.900 CONFIDENCE REGION PLOT ON 400 NM BY 400 NM AREA 



200 
190 
180 
170 
160 
150 
1 40 
130 
120 
110 
100 
90 
80 
70 
60 
50 
40 
30 
20 
10 
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-10 

-20 

-30 

-40 
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-60 

-70 

-80 

-90 

-100 

-110 

-120 

-130 

-140 

-150 

-160 

-170 

-180 

-190 

-200 
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Aft ft *.*• ft \ 
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-200 -160 -120 -80 -40 0 40 80 120 160 200 

CONFIDENCE REGION AREA = 1680 NAUTICAL MILES SQUARED. 



104 



GOLD 


5019.9 


NM 


ON 


AZ 


270.0 


TO 


276.387E 


0.0 . F= 


348.028 


GOLD 


819 .4 


NM 


ON 


AZ 


180.0 


TO 


276. 387F 


-13.737. F = 


15.774 


GOLD 


73.1 


NM 


ON 


AZ 


90.1 


TO 


277.641E 


-13.737. F= 


5.151 


GOLD 


71 . 1 


NM 


ON 


AZ 


180. 0 


TO 


277 . 64 IF 


-14.929. F= 


2 .936 


GRAD 


1 .4 


NM 


ON 


AZ 


267.9 


TO 


277.616E 


-14.930. F= 


2.927 


GRAD 


0.1 


NM 


ON 


AZ 


270.0 


TO 


277.615E 


-14.930. F= 


2.926 



82-23W, 14-56S IS THE LEAST SQUARES SOLUTION. 



OBJECTIVE FUNCTION = 



2.93 DEGREES SQUARED. 



8 DEGRFES OF 


FREEDOM 
















SOLUTION 


SIGNAL 




SITE 








ID INTERCEPT 


AZIMUTH 


A Z I MUT H 




RANGE 


PLOT THETA 


i VALS 


NR SITE NAME 


BETA - 


THETA = OMEGA 


NM 


NORM 


NM 


AZ 


9 HONOLULU 


109 .222 


109.844 - 


0.622 


4943 


204. 5 


36 


294. 5 


17 RENO 


139.373 


138.901 


0.472 


3868 


58.4 


25 


328.4 


6 DETROIT 


179.229 


179. 163 


0.066 


3421 


89.4 


3 


359.4 


14 NASHVILLE 


174.501 


173.573 


0.928 


3060 


84.8 


43 


3 54.8 


3 BOSTON 


192.934 


193.473 - 


0.539 


3479 


280.2 


27 


10.2 


11 MIAMI 


183.252 


184.030 - 


0.778 


2432 


273.6 


30 


3.6 


19 WICHITA 


162 .091 


161 .478 


0.613 


3253 


75.0 


29 


345.0 


8 HELFNA 


148.667 


149.088 - 


0. 421 


4010 


248.4 


23 


338 .4 


MEAN ARC STATISTIC = 59. 


3 DEGREES 












MEAN ACUTE ANGLE OF INTERSECTION 


= 27.8 


DEGREES . 
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GOLD 


9485.0 


NM 


ON 


AZ 


90. 0 


TO 


1 57. 984E 


0.0 . F= 


215.253 


GOLD 


74.4 


NM 


ON 


AZ 


0.0 


TO 


157. 984F 


1.248. F= 


207.475 


GRAD 


1077.6 


NM 


ON 


AZ 


51.3 


TO 


172.232E 


12.394. F= 


6. 162 


GRAD 


13.2 


NM 


ON 


AZ 


15.4 


TO 

TO 


172. 292E 


12.607. F= 


5.864 


GRAD 


0 .0 


NM 


ON 


AZ 


0.0 


172. 292E 


12.607. F= 


5.864 



172-18P, 12-36N IS THE LEAST SQUARES SOLUTION. 

OBJECTIVE FUNCTION = 5.86 DEGREES SQUARED. 

Q DEGREES OF FREEDOM 



SOLUTION SIGNAL SITE 



ID 


INTERCEPT 


AZIMUTH 


AZIMUTH 




RANGE 


PLOT 


THETA 


VALS 


NR 


SI TF NAME 


BETA - 


THETA = 


OMEGA 


NM 


NORM 


NM 


AZ 


9 


HONOLULU 


257.935 


258. 530 


-0.594 


1788 


339.6 


17 


69.6 


10 


LOS ANGELES 


269.320 


268.506 


0.814 


3968 


147. 9 


44 


57.9 


11 


MIAMI 


289.107 


288.596 


0.511 


5995 


150.9 


30 


60.9 


12 


MI DWAY 


213.947 


213. 857 


0. 090 


1098 


120.2 


3 


30.2 


13 


MINNEAPOLIS 


282.050 


282.304 


-0.254 


5064 


315.3 


15 


45.3 


14 


NASHVI LLE 


286. 504 


286.771 


-0.267 


548 0 


322.6 


16 


52.6 


15 


NEW ORLEANS 


284.612 


282.885 


1.727 


5418 


149. 3 


103 


59.3 


16 


NORFOLK 


293.348 


293.874 


-0.526 


5903 


317.3 


31 


47.3 


17 


RENO 


265.894 


266. 982 


-1 . 088 


3907 


322.7 


59 


52 .7 



MFAN ARC STATISTIC = 67.1 DEGREES. 

MEAN ACUTE ANGLE OF INTERSECTION = 13.0 DEGREES. 

0.900 CONFIDENCE REGION PLOT ON 400 NM BY 400 NM AREA 
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CONFIDENCE REGION AREA = 18080 NAUTICAL MILES SQUARED. 
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GOl D 202.5 NM ON 
GOLD 1955.9 NM ON 
GOLD 131'+. 5 NM ON 
GOLD 585.4 NM ON 
GOLD 2030.4 NM ON 
GOLD 56. 1 NM ON 
GRAD 107.2 NM ON 
GRAD 4.2 NM ON 
GRAD 0.1 NM ON 



AZ 90.0 TO 3. 3 73 E 
AZ 0.0 TO 3.373E 
AZ 277.1 TO 337. 319 E 
AZ 0.0 TO 337.319E 
AZ 286.1 TO 290.972E 
AZ 0.0 TO 290. 972 E 
AZ 69.9 TO 293.299E 
AZ 216.4 TO 293.241E 
AZ 215.9 TO 293.240E 



0.0 . F= 

32.769. F = 
32.7 69. F= 
42.536. F= 
42.536. F= 
43.471 . F = 
44.062. F = 
44.006. F = 
44.005. F= 



10268.355 
5296. 391 
3449.257 
2502 .434 
119.602 
51.410 
10.708 
10.553 
10.553 



6 6 —4 6 W ♦ 44- ON IS THE LEAST SQUARES SOLUTION. 
OBJECTIVE FUNCTION = 10.55 DEGREES SQUARED. 
9 DEGREES OF FREEDOM 







SOLUTION 


SIGNAL 




SITE 








ID 


INTERCEPT 


AZI MUTH 


AZIMUTH 




RANGE 


PLOT 


THETA 


> VALS 


NR 


SITE NAME 


BETA - 


THETA = 


OMEGA 


NM 


NORM 


NM 


AZ 


6 


DETROIT 


76.495 


77.395 


-0.900 


721 


178.5 


11 


268. 5 


12 


MI DWAY 


42.761 


41.147 


1.614 


5054 


33.7 


96 


303.7 


18 


SEATTLE 


74.220 


72.903 


1.317 


2294 


24. 5 


48 


294. 5 


5 


DENVER 


68.794 


67 .963 


0.832 


1719 


4. 1 


23 


274.1 


11 


MIAMI 


27. 603 


26.907 


0.696 


1274 


304.7 


15 


214.7 


17 


RENO 


65.683 


66.781 


-1.098 


2357 


193. 2 


41 


283.2 


4 


CORPUS CHRIS 


48.511 


47.035 


1.476 


1765 


335.6 


43 


245.6 


10 


LOS ANGELES 


60.122 


60. 990 


-0.868 


2440 


184.6 


33 


274.6 


16 


NORFOLK 


50.167 


50.470 


-0.303 


534 


146. 7 


2 


236.7 


MEAN ARC STATISTIC = 33. 


6 DEGREES. 











MEAN ACUTE ANGLE OF INTERSECTION = 33.9 DEGREES. 



0.900 CONFIDENCE REGION PLOT ON 400 NM BY 400 NM AREA 
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